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Chapter  1 
Overview 

1.1  Introduction 

This  is  the  final  report  for  Office  of  Naval  Research  Contract  N00014-95-C- 
0124.  Work  under  this  recently  completed  contract  has  focused  on  improving 
the  basic  technology  of  magnetic  resonance  force  microscopy  (MRFM)  and 
working  towards  the  dual  goals  of  dopant  detection  in  silicon  and  single 
electron  spin  detection.  In  addition,  a  number  of  important  advances  were 
achieved,  including  ultrasensitive  force  detection,  magnetic  resonance  charac¬ 
terization  of  dangling  bond  defects  in  Si02,  and  two-  and  three-dimensional 
imaging  of  both  paramagnetic  and  ferromagnetic  resonance  on  the  micron- 
scale. 

Magnetic  resonance  force  microscopy[l,  2,  3]  (MRFM)  is  a  new  scan¬ 
ning  probe  microscope  technique  that  combines  aspects  of  magnetic  res¬ 
onance  imaging  (MRI)  and  atomic  force  microscopy  (AFM).  IBM’s  inter¬ 
est  in  MRFM  is  driven  by  the  possibility  of  achieving  non-invasive,  three- 
dimensional  imaging  with  atomic-resolution  and  elemental  selectivity.  If  this 
goal  can  be  realized,  the  technique  would  have  a  revolutionary  impact  on  the 


1 


2 


CHAPTER  1.  OVERVIEW 


field  of  microscopy  and  have  many  important  applications.  Potential  applica¬ 
tions  include:  1)  determining  the  sub-surface,  three-dimensional  structure  of 
solid  state  materials,  2)  imaging  three-dimensional  distributions  of  dopants 
in  semiconductors  with  angstrom  spatial  resolution,  3)  imaging  defects  and 
trapping  sites  in  semiconductors,  4)  imaging  interactions  between  polymer 
molecules,  5)  determining  the  three-dimensional  atomic  structure  of  macro¬ 
molecules,  such  as  proteins. 

The  long-term  goal  of  single  nuclear  spin  detection  and  imaging  is  an  ex¬ 
tremely  challenging  one  due  to  the  smallness  of  the  forces  that  need  to  be  de¬ 
tected  (e.g.,  attonewtons).  Because  of  the  many  challenges,  IBM  is  pursuing 
a  step-by-step  approach  wherein  new  force  sensing  and  magnetic  resonance 
techniques  axe  being  introduced  in  stages  and  intermediate  results  obtained. 
Microscopes  with  progressively  greater  sensitivity  and  spatial  resolution  are 
constructed  as  new  force  sensing  and  magnetic  resonance  techniques  are  de¬ 
veloped.  Applications  demonstrating  the  efficacy  of  new  techniques  are  being 
developed  in  parallel  with  improvements  in  imaging  capability.  The  recent 
experimental  emphasis  has  been  on  electron  spin  resonance,  instead  of  nuclear 
magnetic  resonance,  because  of  the  larger  magnetic  moment  of  the  electron. 

Work  under  the  ONR  contract  was  aimed  at  extending  the  capability 
of  noicron-scale  MRFM,  while  simultaneously  working  towards  single  spin 
imaging  with  nanometer  resolution.  Specific  accomplishments  include: 

1.  Achievement  of  spin  sensitivity  that  is  100-1000  better  than  conven¬ 
tional  magnetic  resonance  techniques. 

2.  Operation  of  two  low  temperature  MRFM  apparatuses:  one  for  micron- 
scale  studies  and  one  for  single  spin  detection. 

3.  ESR-MRFM  studies  of  a)  phosphorus  dopants  in  silicon,  b)  dangling 
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bonds  in  Si02,  c)  organic  free  radicals  (DPPH),  and  d)  yttrium-iron- 
garnet  ferromagnetic  resonance. 

4.  First  demonstration  of  MRFM  with  a  tip-on-cantilever  configuration. 
True  three-dimensional  (volumetric)  imaging  at  the  micron  scale  was 
demonstrated. 

5.  Demonstration  of  attonewton  force  sensitivity  using  ultrathin  single 
crystal  silicon  cantilevers.  This  is  the  most  sensitive  mechanical  force 
detection  ever  demonstrated. 

6.  Demonstration  of  the  key  ingredients  for  single  spin  MRFM  experi¬ 
ments,  including  a)  operation  of  ultrathin  vertical  cantilevers  within 
10  nm  of  a  sample  surface  at  low  temperature,  b)  development  and 
characterization  of  magnetic  tips  compatible  with  ultrathin  silicon  can¬ 
tilevers;  c)  characterization  of  electron  spins  in  a  test  sample  consisting 
of  dangling  bond  defects  in  Si02  (E'  centers). 

7.  Demonstration  of  quantitative  magnetic  field  mapping  using  force-detected 
ferromagnetic  resonance. 


1.2  Organization  of  this  report 

This  report  is  organized  into  chapters  that  describe  various  aspects  of  MRFM 
covered  under  the  ONR  contract.  Portions  of  the  report,  especially  chapters 
2-5,  are  adapted  from  the  Ph.D.  thesis  of  Koichi  Wago,  who  was  partially 
supported  by  the  contract. 

Chapter  2  overviews  basic  principles  of  low  temperature  MRFM  and  ex¬ 
perimental  procedures.  Force  resolution  achievable  at  low  temperature  is  also 
discussed  in  Chapter  2. 


4 


CHAPTER  1.  OVERVIEW 


Chapter  3  describes  the  detection  of  phosphorus  dopants  in  silicon  by 
force-detected  electron  paramagnetic  resonance.  We  also  discuss  experiments 
on  force-detected  EPR  spectroscopy,  including  the  observation  of  the  hyper- 
fine  splitting  due  to  nuclei. 

Chapter  4  describes  EPR  measurements  of  E'  centers  (dangling  bond 
defects)  in  Si02,  including  spectroscopy,  spin  nutation  and  spin  echoes.  Re¬ 
laxation  times  of  electron  spins  while  undergoing  cyclic  adiabatic  inversion 
are  presented  and  discussed  in  terms  of  the  utility  for  future  single-spin  de¬ 
tection  experiments. 

Chapter  5  describes  an  improved  MRFM  configuration  incorporating  a 
magnetic  tip  mounted  directly  on  the  cantilever.  The  improved  microscope 
was  used  to  compare  and  contrast  two-  and  three-dimensional  imaging  char¬ 
acteristics  of  EPR  and  ferromagnetic  resonance  (FMR).  Quantitative  mag¬ 
netic  field  mapping  on  the  micrometer  scale  is  also  discussed. 

Chapter  6  discusses  the  ongoing  effort  aimed  at  detecting  individual  elec¬ 
tron  spins.  It  describes  the  basic  experimental  approach  and  the  status  of 
key  experimental  sub-systems,  including  fabrication  of  ultrasensitive  silicon 
cantilevers  capable  of  detecting  attonewton  forces. 


1.3  Bibliography  of  papers  published  acknowl¬ 
edging  ONR  support 
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Chapter  2 

Basic  Principles 

2.1  Why  low-temperature  MRFM? 

Improving  the  detectability  of  small  forces  is  a  key  aspect  of  MRFM  research 
since  the  force  resolution,  in  combination  with  the  magnitude  of  the  tip 
field  gradient,  sets  the  ultimate  spin  sensitivity.  The  force  resolution  of  a 
cantilever-based  force  sensor  is  fundamentally  limited  by  thermomechanical 
noise.  If  a  cantilever  with  a  spring  constant  k  is  in  thermal  equilibrium  at 
temperature  T,  the  equipartition  theorem  predicts 

ifc(x^)  =  (2.1) 

where  (x^)  is  the  mean-square  vibration  amplitude  of  the  cantilever  and  ks 
is  the  Boltzmann  constant.  If  we  assume  that  the  vibration  noise  arises  from 
a  force  noise  with  a  fiat  power  spectrum  Sp.,  then  the  resulting  cantilever 
vibration  amplitude  is 

{x^)  =  (2.2) 
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where 


GH  = 


—  a;2)  +  iujcO^IQ 


(2.3) 


is  the  transfer  function  of  the  cantilever  with  a  spring  constant  k,  resonant 
frequency  Wc,  quality  factor  Q.  Integration  of  Eq.  (2.2)  for  Q  ^  1  gives 


(x^)  =  Sf 


<^cQ 
4fc2  ' 


Prom  Eqs.  (2.1)  and  (2.4),  we  obtain 


Sf  = 


AkkgT 

^cQ 


(2.4) 


(2.5) 


and  the  minimum  detectable  force  in  a  detection  bandwidth  Ai/  is  given 
by[4,  5,  6] 

(2.6) 

It  is  clear  from  Eq.  (2.6)  that  low-temperature  operation  should  increase 
the  sensitivity  of  MRFM.  An  additional  advantage  of  low-temperature  opera^- 
tion  is  that  the  magnetic  moment  due  to  spin  polarization,  which  is  inversely 
proportional  to  temperature,  will  also  increase,  resulting  in  a  larger  force  sig¬ 
nal.  Equation  (2.6)  also  suggests  the  strategy  on  choosing  cantilevers:  use  a 
soft  (small  k)  cantilever  with  high  Q.  Operating  at  low  temperature  will  ben¬ 
efit  in  this  sense  as  well,  since  Q  generally  improves  at  lower  temperature[7] . 

This  chapter  and  chapters  3-5  concentrate  on  low  temperature  MRFM 
experiments  using  commercially  available  cantilevers  capable  of  force  detec¬ 
tion  in  the  femtonewton  range.  Chapter  6  shows  how  the  force  resolution 
can  be  extended  to  the  attonewton  range  using  ultrathin  custom-fabricated 
cantilevers. 
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2.2  Low-temperature  MRFM  configuration 

Figure  2.1  shows  the  schematic  of  an  experimental  apparatus  used  for  micron- 
scale  MRFM  experiments.  The  MRFM  apparatus  was  suspended  in  a  vac¬ 
uum  can  in  a  liquid  helium  dewer[8]  with  a  superconducting  magnet  to  gen¬ 
erate  a  polarizing  magnetic  field.  The  apparatus  was  cooled  via  a  copper 
braid  connected  to  the  bottom  of  the  can  to  provide  thermal  conduction 
while  maintaining  adequate  vibration  isolation. 

To  ensure  that  the  resonant  volume  is  large  enough  to  generate  a  mea¬ 
surable  signal,  a  millimeter-size  magnetic  tip  was  placed  external  to  the  can¬ 
tilever  and  the  sample  was  directly  glued  onto  the  cantilever.  For  example, 
an  iron  tip  of  1.5  mm  diameter  and  3  mm  long  was  placed  roughly  700  jum 
away  from  the  cantilever  for  the  NMR  experiments  in  Chapter. 

The  rf  field  Bi  was  applied  to  the  sample  using  a  nearby  coil  to  induce 
magnetic  resonance.  The  polarizing  magnetic  field  Bq  was  a  combination  of 
the  field  from  the  superconducting  magnet  and  the  field  from  the  magnetic 
tip.  The  oscillatory  magnetic  force  resulting  from  magnetic  resonance  excites 
the  vibration  of  the  cantilever.  The  angstrom-scale  motion  of  the  cantilever 
was  measured  with  a  fiber-optic  interferometer  [9]  and  was  synchronously 
detected  with  a  lock-in  amplifier. 

It  might  seem  that,  since  the  magnetic  force  is  proportional  to  the  field 
gradient,  using  a  smaller  tip  placed  closer  to  the  sample  should  increase  the 
signal  due  to  the  stronger  field  gradient.  In  general,  scaling  down  dimensions 
of  the  tip  size  and  tip-sample  distance  increases  the  field  gradient  in  inverse 
proportion  to  the  scale  factor.  This  has  the  beneficial  effect  of  increasing  the 
force  per  spin  (or  per  unit  volume).  However,  the  resonant  volume  decreases 
faster  as  the  dimensions  axe  scaled  down,  because  not  only  the  thickness  of 
the  resonant  slice  scales  down  in  proportion  to  the  scale  factor,  but  the  area 
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fiber  optic 
interferometer 


superconducting 

magnet 


I1mffrn2.cdr 


Figure  2.1:  Low  temperature  MRFM  setup. 
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covered  by  the  resonant  slice  decreases  as  well.  As  a  result,  the  net  force, 
(force  per  volume)  x  (resonant  volume),  decreases  as  the  dimensions  are 
scaled  down.^ 

In  the  configuration  shown  in  Fig.  2.1,  where  a  micron-size  sample  and  a 

millimeter-size  magnetic  tip  are  used,  the  resonant  volume  is  given  by 

A  n 

where  Ax  and  Ay  are  the  sample  lateral  dimensions  and  AB  is  the  magnetic 
resonance  linewidth.^  The  force  per  volume  is  given  by 

F  =  (2.8) 

where  N  is  the  number  of  spin  per  volume.  The  net  force  is  then  given  by 

Fnct  =  '^^°Aa;AyAB.  (2.9) 

Substituting  typical  numbers^  N  =  4.9  x  10^^  cm~^,  /i  =  1.33  x  10“^®  J/T, 
Bo  -  2.5  T,  Ax  =  Ay  -  10  /xm,  ks  =  1.38  x  10“^^  J/K,  T  =  5  K,  and 
AB  =  2  G,  we  obtain  F„et  6  x  10“^^  N,  which  should  be  detectable  as 
discussed  in  the  next  Section. 


2.3  Force  resolution 

For  experiments  described  in  Chapters  2-6,  we  used  commercial  single-crystal- 
silicon  cantilevers[10]  with  dimensions  of  approximately  470  /ixm  X  45  ^m  X 

^This  argument  does  not  apply  to  the  case  of  single-spin  MRFM,  for  which  the  gradient 

may  be  increased  without  concern  for  the  size  of  the  resonant  volume. 

^The  sample  is  assumed  to  be  thicker  than  the  resonant  slice  such  that  Az  > 

AR/|ViBtip|. 

^The  values  of  N  and  are  for  in  1%  Nd-doped  CaF2.  Linewidth  AB  is  generally 
determined  by  the  natural  linewidth  and  the  modulation  amplitude  used  to  generate  the 
oscillatory  force  as  discussed  later. 
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1.5  //m,  a  spring  constant  k  ~  0.07  N/m,  and  a  resonance  frequency  fc  —  9.8  kHz. 
Figure  2.3  shows  the  optical  micrograph  of  a  typical  cantilever  with  sample 
loaded.^ 

The  single-crystal-silicon  cantilevers  showed  very  high  Q  in  vacuum 
10“®  torr)  and  at  low  temperature.  Q  was  determined  by  measuring  the 
decay  time  of  the  vibration  amplitude  after  the  vibration  was  driven  by  a 
piezoelectric  actuator  at  the  base  of  the  cantilever.  Figure  2.3  shows  a  typi¬ 
cal  Q  measurement  for  a  cantilever  with  sample  loaded  (resonant  frequency 
fc  =  1.6  kHz).  Decay  time  becomes  longer  as  the  temperature  is  decreased, 
yielding  Q  ~  500  000  at  6  K,  improved  from  Q  ~  20  000  at  room  temperature. 

Low  temperature  force  resolution  was  characterized  using  a  bare  silicon 
cantilever  {fc  =  9.8  kHz).  Figure  2.4  shows  the  thermal  vibration  ampli¬ 
tude  of  the  cantilever  measured  synchronously  at  the  cantilever  resonant 
frequency.  The  vibration  amplitude  was  1.1  A-rms  at  room  temperature 
and  decreased  to  0.26  A-rms  at  6  K.  Prom  the  spring  constant,  the  Q 
at  each  temperature,®  and  the  position  of  the  optical  fiber,  we  estimate® 

■^The  samples  were  glued  onto  cantilevers  with  a  UV-curing  epoxy  (Loctite).  A  pipette 
with  a  sharp  pulled  end  (prepared  using  lighter  flame)  and  a  micromanipulator  were  used 
under  a  microscope  to  mount  the  sample  on  the  cantilever. 

~  40  000  at  room  temperature  and  Q  ~  200  000  at  6K  for  this  cantilever. 

®The  procedure  we  used  to  determine  Sp^  is  as  follows.  The  root-mean-square  can¬ 
tilever  vibration  amplitude,  A,  was  measured  at  both  300  K  and  6  K  using  a  lock-in  am¬ 
plifier  with  a  short  time  constant  (rx,  =  100  ms).  The  short  time  constant  (tjt,  •C  Q/wc) 
ensured  that  the  entire  bandwidth  of  the  cantilever  noise  was  detected.  A  correction  fac¬ 
tor,  c,  was  applied  to  account  for  the  fact  that  the  optical  fiber  was  not  measuring  the 
motion  of  the  cantilever  end,  but  rather  at  a  point  approximately  40%  from  the  end. 

The  correction  factor  was  determined  using  room-temperature  thermal  noise  such  that 
A(cA3oojr)^/2  =  kBT/2  for  T  =  300  K.  The  determined  value,  c  =  2.3,  was  consistent 
with  that  estimated  from  the  mode  shape,  4>{^)  —  [cosh  k®— cos  /c®— ©"(sinh  /c®— sin  k®)]/2, 
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Single-crystal  Si  cantilever 


Sample 


Figure  2.2:  An  optical  micrograph  of  a  typical  cantilever  and  sample. 


Vibration  amplitude  (a.  u.) 
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Figure  2.3:  Cantilever  vibration  decay.  The  motion  of  the  cantilever  was 
monitored  by  the  fiber-optic  interferometer  and  the  amplitude  of  the  vibra¬ 
tion  at  the  resonant  frequency  was  measured  after  the  vibration  was  driven. 


14 


CHAPTER  2.  BASIC  PRINCIPLES 


Sp!^  =  8  X  10-^^  N/\^  at  6  K.  This  is  a  factor  of  8  improvement  com¬ 
pared  to  room  temperature.  The  measured  force  resolution  is  about  two 
times  larger  than  the  theoretical  value,  4.3  X  10"^^  N / -s/Hz,  calculated  from 
Eq.  (2.5).  The  discrepancy  is  probably  due  to  some  environmental  vibration 
noise. 


2.4  Active  damping  feedback 


In  addition  to  the  improved  force  resolution,  the  high  Q  of  the  silicon  can¬ 
tilever  also  means  that  it  takes  a  long  time  for  the  vibration  amplitude  to 
reach  the  steady  state  value.  In  order  to  speed  up  the  cantilever  response,  an 
electromechanical  feedback  technique  was  often  used[ll,  12,  13].  A  damp¬ 
ing  feedback  signed,  90°  phase-shifted  from  the  measured  cantilever  vibration 
signal  x(t),  was  applied  to  the  piezoelectric  actuator  at  the  base  of  the  can¬ 
tilever. 

The  feedback  force  is  proportional  to  x{t)  and  can  be  written  as  Ffeedbaefc(tt’) 
—iau}x[u})  in  the  Fourier  domain,  where  cc  is  a  constant.  The  cantilever  vi¬ 
bration  amplitude  x{oj)  driven  by  a  force  F{lo)  is  then  described  as 

x{u})  =  G{u})  [F(w)  -I-  Ffeedback(w)]  =  G{uj)  [F{w)  -  iaujx{u})] ,  (2.10) 


where 


G(a;)  = 


wHk 


(a;2  -  uA)  -f  iujciolQ 
is  the  transfer  function  of  the  cantilever  as  in  Eq.  (2.3). 
Eq.  (2.10)  as 

x{oj)  =  G^s{u})F{(jj), 


(2.11) 

We  can  rewrite 

(2.12) 


where  k  1.875/jL  and  a  0.734.  Spectral  density  of  the  force  noise  at  6  K  was  estimated 
according  to  5^^  =  UcA^kIQ  •  (units  of  N/v1Sz).  The  quantity  ^Q/oJc  is  the 

inverse  of  the  equivalent  noise  bandwidth  of  the  cantilever  resonance. 
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(b) 


Figure  2.4:  The  vibration  noise  of  the  cantilever.  The  motion  of  the  cantilever 
was  monitored  by  the  fiber-optic  interferometer  and  the  amplitude  of  the 
vibration  at  the  resonance  frequency  was  measured  at  (a)  300  K  and  (b)  6 
K.  The  measurement  bandwidth  was  determined  by  the  natural  bandwidth 
of  the  cantilever,  u>cl4Q  =  0.4  Hz  for  T  =  300  K  and  0.08  Hz  for  6  K. 
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where 

r  (u)  - 

with  an  effective  Q 


(2.13) 


(2.14) 


As  a  result,  the  1/e  response  time  of  the  cantilever  is  reduced  to  Tes  = 
2Qeff/a;c  from  r  =  2Qlujc- 

This  reduction  of  effective  Q  and  r  occurs  without  degrading  the  force 
resolution  for  the  same  detection  bandwidth,  since  the  active  damping  does 
not  affect  the  spectral  density  of  the  noise  force  Sp.  Sp  is  determined  by  the 
coupling  to  the  thermal  bath  and  given  by  Sp  =  AkkpT/wcQ  as  in  Eq.  (2.5), 
where  Q  is  the  intrinsic  Q  of  the  cantilever. 

The  advantage  of  using  feedback  to  speed  up  the  cantilever  response 
rather  than  using  a  cantilever  with  low  intrinsic  Q  can  be  understood  if 
we  consider  the  energy  of  the  actively  damped  cantilever.  The  mean-square 
vibration  amplitude  due  to  thermal  noise  is  given  by 


(l=)  =  ^  SF|G.*(a,)|^^ 


<3eff 


rspiGiu)]-^ 

Jo 


duj 


(2.15) 


and  is  therefore  reduced  by  a  factor  of  QctijQ  compared  to  the  case  without 
feedback  given  by  Eq.  (2.2).  The  energy  of  the  cantilever  is 

\k{x'‘)  =  ^  .  \hBT,  (2.16) 

and  no  longer  in  thermal  equilibrium  with  its  environment. 
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2.5  Techniques  for  generating  oscillatory  mag¬ 
netic  forces 

The  maximum  vibrational  response  of  the  cantilever  is  achieved  when  the 
sample  magnetization,  and  hence  the  resulting  magnetic  force,  is  cyclically 
modulated  at  the  mechanical  resonant  frequency  of  the  cantilever  Wg.  Typ¬ 
ically  iVc  is  in  the  kilohertz  range,  much  slower  than  the  Larmor  frequency. 
Thus,  at  least  for  our  present  cantilevers,  the  precession  of  the  magnetization 
cannot  be  directly  used  as  the  source  of  oscillatory  forces  matched  to  a>c.  In¬ 
stead,  two  techniques,  described  below,  have  been  used  to  modulate  the  lon¬ 
gitudinal  component  of  the  magnetization  M^,  causing  Fz  =  Mz{dBtipldz) 
to  oscillate  at  a>c. 


2.5.1  Cyclic  saturation 


In  general,  the  behavior  of  the  magnetization  is  described  phenomenologically 
by  the  Bloch  equation  given  by[14] 


dM 

dt 


=  M  X  7B  — 


A4i  -F  Myj  {Mz  -  Mo)k 
T2  Ti 


(2.17) 


where  i,  j,  and  k  are  unit  vectors  in  the  x,  y,  and  z  directions,  respectively, 
and  Mo  is  the  equilibrium  magnetization,  Ti  and  T2  are  spin-lattice  and  spin- 
spin  relaxation  times,  respectively.  In  the  presence  of  the  static  field  Bq  in  the 
z  direction  and  the  rf  field  Bi  at  frequency  a;  in  the  x  direction,  it  is  customary 
to  introduce  the  rotating  frame  at  frequency  w,  where  Bi  is  stationary.^  The 
steady  state  solution  of  the  Bloch  equation  (i.e.,  dM-ldt  =  0)  in  the  rotating 


is  decomposed  into  two  circularly  polarized  components  rotating  in  opposite  di¬ 
rections  in  the  xy  plane  and  we  use  the  approximation  to  ignore  the  counter  rotating 
component. 
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frame  leads  to 


M,  =  Mo 


T^BlTiT2 


(2.18) 


i  +  i'rBo-ujyTi  +  'r^BlTiT2\' 

For  the  spins  that  are  excited  at  resonance  (w  =  'fBo),  Mz  is  suppressed  to 
the  value  Mq  —  SMz,  where 

=  Twmr.- 

Modulating  either  the  amplitude  Bi,  frequency  oj  of  the  rf  field,  or  the 
polarizing  field  strength  Bq  at  frequency  lOc  cyclically  suppresses  Mz,  causing 
the  magnetic  force  Fz  to  oscillate  at  lOc-  An  assumption  made  here  is  that  the 
magnetization  reaches  the  steady  state  value  sufficiently  fast  (Ti,  T2  1/wc)? 
which  is  typically  the  case  for  electron  spin  systems.  This  technique  was  used 
for  the  first  EPR-MRFM  experiments  [2,  15,  16]. 


2.5.2  Cyclic  adiabatic  inversion 

For  spins  with  long  relaxation  times,  the  cyclic  saturation  technique  is  not 
effective.  An  alternative  method  is  the  cyclic  adiabatic  inversion  technique 
used  in  NMR-MRFM[11,  17]. 

In  the  rotating  frame,  Bq  is  transformed  to  Bq  —  ujIj  as  shown  in  Fig. 
2.5[18].  The  magnetization  processes  about  the  effective  field 

Beff  =  BiX+  (2.20) 

where  x  and  z  are  unit  vectors  in  the  rotating  frame.  Beff  can  be  manipulated 
by  modulating  the  frequency  oj.  If  the  direction  of  Beff  is  changed  starting 
from  being  parallel  to  the  magnetization,  then  the  magnetization  follows  the 
direction  of  Beff  as  long  as  the  rate  of  change  is  “slow”  enough  such  that 

dBeff 

dt 


« iBl 


(2.21) 
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laboratory  frame 


rotating  frame 
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Figure  2.5:  The  laboratory  frame  and  rotating  frame. 
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This  condition  is  described  by  the  term  “adiabatic”  condition[18]. 

For  example,  modulating  the  frequency  of  the  rf  field  uj  such  that 

u}{t)  =  (jjQ  —  Qcosojct,  (2.22) 

where  (jJq  =  7B0,  causes  Bes  to  swing  up  and  down  at  ujc-  If  the  adiabatic 
condition  is  satisfied,  the  magnetization  is  cyclically  inverted  following  Beff, 
resulting  in  the  oscillatory  magnetic  force  F^.  In  terms  of  frequencies,  the 
adiabatic  condition  (2.21)  becomes 

<  (2.23) 

or 

fiwc  <  7"^1-  (2.24) 


* 


Chapter  3 

MRFM  Detection  of  Dopants 
in  Silicon 


One  potential  technological  application  of  MRFM  is  the  imaging  of  dopants 
in  semiconductors.  This  chapter  describes  basic  experiments  that  demon¬ 
strate  detection  and  spectroscopy  of  phosphorus  dopants  in  silicon  using 
force-detected  electron  paramagnetic  resonance  (EPR).^  Force-detected  EPR 
spectra  were  obtained  using  an  amplitude  or  frequency  modulated  microwave 
field  to  cyclically  saturate  the  spin  magnetization.  For  a  sample  containing 
4  X  10^®  phosphorus  atoms/cm^,  a  single  strong  EPR  line  was  observed.  For  a 
sample  containing  8  x  10^®  phosphorus  atoms/cm^,  a  pair  of  lines  split  by  the 
42-gauss  ^^P  hyperfine  interaction  was  observed.  This  result  demonstrates 
the  possibility  of  using  MRFM  techniques  for  spectroscopic  purposes. 


portion  of  the  work  presented  in  this  chapter  was  performed  after  the  submission 
of  the  ONR  contract  proposal  in  1995,  but  before  the  official  start  of  the  contract.  We 
include  the  results  here  since  they  are  directly  relevant  to  the  goals  of  the  contract 
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Fiber-optic 

interferometer 


Figure  3.1:  Schematic  setup  of  the  experiment 

3.1  Experiment 

The  experimental  apparatus  is  shown  schematically  in  Fig.  3.1. 

Two  samples  with  different  phosphorus  concentrations  were  studied.  The 
“high-concentration”  sample  had  a  phosphorus  concentration  No  =  4  x 
10^®  cm“®,  and  had  a  size  of  roughly  25/xm  x  25/im  x  5/im.  The  “low- 
concentration”  sample  had  =  8  x  10^®  cm"®  and  was  slightly  larger  in  size. 
The  samples  were  glued  onto  commercial  single-crystal-silicon  cantilevers [10] 
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with  spring  constant  k  ~  0.07  N/m.  The  resulting  mechanical  resonant 
frequencies  fc  were  7.6  kHz  and  5.0  kHz  for  the  high  and  low  concentra¬ 
tion  samples,  respectively.  The  vibration  amplitude  of  the  cantilever  was 
monitored  with  a  fiber-optic  interferometer[9]  using  an  optical  power  of  less 
than  10  /iW  at  the  cantilever.  At  higher  optical  power,  we  found  that  the 
cantilever  Q  was  affected  by  an  interaction  between  the  laser  light  and  the 
cantilever  and,  depending  on  the  condition  of  optical  interference,  Q  was 
either  suppressed  or  enhanced. 

EPR  was  excited  by  a  microwave  field  Bi  at  a  frequency  a;/27r  =  12.6  GHz 
using  an  end-coupled  half-wavelength  microstripline  resonator [19,  20].  The 
resonator  consisted  of  a  4  mm  long  section  of  50-ohm-impedance  microstrip 
transmission  line  fabricated  on  a  0.5  mm  thick  sapphire  substrate.  With  an 
input  power  of  24  mW,  the  resonator  produced  0.6  G  at  the  sample  position 
200  fj,ra  above  the  center  of  the  resonator.^  A  NdFeB  permanent  magnet,  3 
mm  in  both  diameter  and  length,  was  located  800  /xm  away  from  the  sample, 
and  produced  a  field  gradient  G  =  2  GJixm.  This  served  to  generate  a 
magnetic  force  on  the  sample,  as  well  as  defining  a  “resonant  slice”  within 
which  the  magnetic  resonance  condition,  u)  =  'yBo,  is  satisfied.  The  total 
field  Bq  was  a  combination  of  the  fields  from  the  permanent  magnet  (~  2000 
G)  and  a  superconducting  magnet.  A  control  experiment  was  performed 
to  verify  that  antimony  dopants  in  the  cantilever  contribute  no  EPR  force 
signal.  Presumably  this  is  because  of  the  very  short  relaxation  times  of 
antimony-doped  silicon[21]. 


^The  strength  of  Bi  was  calibrated  by  spin  nutation  experiments  described  in  Chapter 
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3.2  Force  detection  of  electron  paramagnetic 
resonance 

Cyclic  saturation,  described  in  Section  2.5.1,  was  used  to  induce  an  oscil¬ 
latory  component  in  the  longitudinal  magnetization  M^.  Either  amplitude 
modulation  of  Bi  (AM  cyclic  saturation)  or  frequency  modulation  of  Bi  (FM 
cyclic  saturation)  was  used,  as  described  below. 


3.2.1  AM  cyclic  saturation 


For  the  case  of  AM  cyclic  saturation,  Bi  was  modulated  at  frequency  fc  so 
that  Bi  was  turned  on  and  off  with  a  50%  duty  cycle.  As  a  result,  oscil¬ 
lates  between  the  values  Mo  and  Mq  —  SM^  (see  Eq.  (2.19))  with  frequency 
/c,  provided  that  Ti  is  short  enough  so  that  Mz  reaches  its  steady  state  value 
quickly  compared  with  the  time  scale  of  1/ fc-  The  oscillating  magnetization 
has  the  first-harmonic  Fourier  component  Mi  =  {2I'k)5Mz,  which  produces 
an  oscillating  magnetic  force  MiG  in  the  presence  of  the  field  gradient  G. 
Spins  that  axe  slightly  off  resonance  also  contribute  to  the  oscillatory  force 
but  less  effectively  since  the  suppression  of  Mz  is  not  complete.  The  total 
force  is  obtained  by  integrating  these  contributions  over  the  sample  volume 
taking  into  account  that  Bq  is  a  function  of  position  z.  Therefore 


Mo't^BlTiT2 

1  +  bBo{z)  -  uj^T^  +  'r^BlTiT^ 


GA{z)  dz, 


(3.1) 


where  A{z)  is  the  cross-sectional  area  of  the  sample  at  z.  The  dominant  con¬ 
tribution  to  the  total  force  originates  from  the  spins  in  the  resonant  slice,  the 
thickness  of  which  is  on  the  order  of  AB/ G,  where  AB  =  2\J\  -\-  '^^BiTiT^I  (7^2) 
is  the  full-width-at-half-maximum  of  the  resonance  given  by  Eq.  (2.18).  If 
the  strength  of  the  polarizing  field  Bq  is  scanned  by  changing  the  external 
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field  from  the  superconducting  magnet,  the  oscillatory  magnetic  force  F  is 
expected  to  be  at  its  maximum  when  the  resonant  slice  falls  roughly  in  the 
middle  of  the  sample. 

Figure  3.2(a)  shows  the  force-detected  EPR  spectrum  for  the  high-concentration 
sample  at  T  =  5.5  K.  Bi  amplitude  modulation  wa^  applied  and  Bq  was  swept 
while  measuring  the  cantilever  vibration  amplitude,  A  single  peak  was  ob¬ 
served  with  a  maximum  force  of  about  N.  The  observed  linewidth  of  6  G 

is  not  the  intrinsic  linewidth  AB,  but  rather  includes  additional  broadening 
due  to  the  extended  sample  thickness  in  the  presence  of  the  field  gradient. 


3.2.2  FM  cyclic  saturation 

For  the  case  of  FM  cyclic  saturation,  the  microwave  frequency  was  mod¬ 
ulated  as  a;(t)  =  wq  +  fisin(27r/cf),  generating  a  time  varying  Mz(t)  that 
contains  the  first-harmonic  Fourier  component  Mi.  For  small  values  of  fi. 
Ml  is  given  by  DdM^Idu).  Mi  is  zero  exactly  on  resonance  {u>o  =  'fBo)  and 
has  positive  and  negative  peaks  at  Bq  =  a;o/7  d=  AB  1 2-\/Z  when  Bq  is  swept 
through  the  resonance.  For  an  extended  sample,  the  bipolar  nature  of  the 
resonance  line  leads  to  cancellation  of  the  net  force  when  the  resonant  slice 
is  in  the  center  of  the  sample.  The  net  force  is  maximum  when  the  reso¬ 
nant  slice  is  positioned  at  either  edge  of  the  sample.  Figure  3.2(b)  shows 
the  force-detected  EPR  spectrum  measured  with  FM  cyclic  saturation  for 
the  high-concentration  sample.  The  peak  deviation  J2/27r  was  10  MHz  and 
T  =  6.5  K.  The  spectrum  shows  positive  and  negative  peaks,  resembling  the 
derivative  of  Eq.  (2.18)  as  expected.  The  observed  linewidth  was  6  G.  There 
are  three  factors  contributing  to  this  linewidth:  (1)  the  intrinsic  linewidth 
of  the  resonance,  (2)  broadening  due  to  the  frequency  modulation,  and  (3) 
broadening  due  to  the  sample  thickness  in  the  presence  of  the  field  gradient. 
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Figure  3.2:  Force-detected  EPR  spectrum  of  the  high-concentration  sample, 
(a)  AM  cyclic  saturation  with  T  =  5.5  K.  (b)  FM  cyclic  saturation  with 
T  =  6.5K  and  fi/27r  =  10  MHz. 
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The  difference  between  the  positive  and  the  negative  force  peaJc  is  plotted 
as  a  function  of  temperature  in  Fig.  3.3.  The  force-detected  EPR  signal 
increased  as  the  temperature  decreased  and  was  approximately  proportional 
to  Due  to  the  nature  of  degenerate  conduction-electron  paramagnetism 

(Pauli  paramagnetism),  the  temperature  dependence  of  the  magnetization 
Mo  can  be  neglected  in  this  temperature  range[22].  We  believe,  therefore, 
that  the  temperature  dependence  of  the  force-detected  EPR  signal  is  mainly 
due  to  the  changing  relaxation  times,  Ti  axid  T^.  Assuming  that  the  resonant 
slice  is  positioned  at  the  edge  of  the  sample  (z  =  0),  and  that  the  sample 
is  much  thicker  than  the  resonant-slice  thickness  and  has  a  constant  cross- 
sectional  area  A,  we  obtain  the  peak  EPR  force  as 


=  r 

Jo 


duj 


l^BlTiT2 


1  -h  {7[^o(0)  -f-  Gz]  -  u^Ti  + 


GAdz 


=7Bo(0) 


DMoA'yBlTiT2 
l+'r^B^TiT2  ■ 


(3.2) 


For  Bi  values  that  are  not  too  large,  F  is  expected  to  increase  as  the  re¬ 
laxation  times  become  longer  at  lower  temperatures  [23],  consistent  with  the 
result  in  Fig.  3.3. 


3.2.3  Observation  of  hyperfine  splitting 

For  the  low-dopant-concentration  sample,  additional  features  were  observed 
in  the  spectrum.  Figure  3.4  shows  the  force-detected  EPR  spectrum  for  FM 
cyclic  saturation  at  T  =  10  K.  The  peak  deviation  Q,/2ir  was  10  MHz  and  the 
signal  was  averaged  over  eight  Bq  sweeps  to  increase  the  signal-to-noise  ratio. 
The  spectrum  shows  two  separate  bipolar  EPR  lines  spaced  by  42  G.  The 
origin  of  the  two  lines  is  the  hyperfine  interaction  between  the  electron  spins 
and  phosphorus  nuclear  spins  [24].  It  is  well  known  that  at  low  temperature 
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Figure  3.3:  Temperature  dependence  of  the  force-detected  EPR  signal  of  the 

high-concentration  sample.  For  reference,  the  dotted  line  shows  a  slope  of 
y-1.5 
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Figure  3.4;  Force-detected  EPR  spectrum  of  the  low-concentration  sample 
using  FM  cyclic  saturation  for  T  =  10  K  and  Q/2w  =  10  MHz.  The  42-gauss 
spacing  between  the  two  bipolar  features  is  due  to  hyperfine  splitting. 
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and  at  low  doping  concentration,  the  unpaired  electron  spin  from  the  donor 
atom  is  localized  on  the  donor  atom,  allowing  the  hyperfine  splitting  to  be 
revealed[25].  (In  contrast,  for  the  high-concentration  sample,  the  observed 
EPR  corresponds  to  conduction  electrons[22].) 

The  peak-to-peak  force  is  plotted  as  a  function  of  temperature  in  Fig. 
3.5.  For  temperatures  above  10  K,  the  EPR  force  increased  rapidly  as  the 
temperature  decreased.  It  showed  stronger  temperature  dependence,  ~ 
than  for  the  high-concentration  sample.  This  is  likely  due  to  stronger  tem¬ 
perature  dependence  of  both  Mq  and  the  relaxation  times.  Unlike  the  high- 
concentration  sample,  the  paramagnetism  of  electrons  localized  on  donor 
atoms  obeys  the  Curie-law  temperature  dependence,  Mq  ~  T~^.  It  is  also 
known  that  Ti  increases  as  ~  T"",  where  n  is  in  the  rage  of  7  to  9,  at 
low  temperature[26,  27].  The  strong  temperature  dependence  of  Mq  and  Ti 
results  in  the  rapid  increase  of  the  force  as  expected  from  Eq.  (3.2).  For 
temperatures  below  10  K,  the  EPR  force  decreased  as  the  temperature  de¬ 
creased.  This  cannot  be  explained  simply  from  Eq.  (3.2).  However,  in  this 
temperature  range,  Ti  becomes  comparable  to  the  time  scale  of  l//c[26,  27], 
which  prevents  the  magnetization  from  reaching  the  steady  state  value  in  Eq. 
(2.18)  when  it  is  modulated  at  fc.  This  results  in  a  smaller  EPR  force  than 
expected  by  Eq.  (3.2),  suppressing  the  signal  at  lower  temperatures. 

3.2.4  Interrupted  cyclic  adiabatic  inversion 

When  Ti  is  long  compared  with  1/fc,  cyclic  saturation  is  not  an  effective 
way  to  generate  the  cantilever  vibration.  In  this  situation,  an  alternative 
way  to  drive  the  cantilever  is  using  interrupted  cyclic  adiabatic  inversion[28]. 
This  is  performed  by  applying  frequency  and  amplitude  modulation  of  Bi 
as  shown  in  Fig.  3.6.  During  the  inversion  cycle,  the  magnetization  is  in- 
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Figure  3.5:  Temperature  dependence  of  the  force-detected  EPR  signal  for  the 
low-concentration  sample.  For  reference,  the  dotted  line  indicates  a  slope  of 
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Figure  3.6:  Amplitude  and  frequency  modulation  pattern  used  for  inter¬ 
rupted  cyclic  adiabatic  inversion. 

verted  down  and  then  up  by  manipulating  the  direction  of  the  effective  field 
Befr(t)  =  [Bo  — a;(t)/7]z-fBiX  in  the  rotating  frame[28]^  (z  and  x  are  the  unit 
vectors  in  the  rotating  frame).  This  generates  an  oscillatory  magnetic  force 
on  the  cantilever.  Between  inversion  cycles,  a  repolarization  period  lasting 
two  cantilever  cycles  is  inserted  to  ensure  that  the  magnetization  does  not 
die  out  over  many  cycles  of  inversion.  This  method  of  manipulating  the  mag- 


®See  Section  2.5.2 
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netization  works  only  if  Ti  is  comparable  to  or  longer  than  1/ fc]  if  Ti  is  too 
short,  FM  cyclic  saturation  results.  The  two  cases  should  be  distinguished 
by  different  shapes  of  the  EPR  force  spectrum,  since  the  cyclic  adiabatic 
inversion  has  a  unipolar  spectral  response  with  the  maximum  force  obtained 
when  the  resonant  slice  is  in  the  middle  of  the  sample.  Figure  3.7  shows 
the  EPR  force  spectrum  for  the  low-concentration  sample  obtained  with  this 
technique  at  T  =  6.5  K  with  Q,/2tt  —  10  MHz.  The  spectrum  was  averaged 
over  eight  sweeps.  It  shows  two  peaks  separated  by  42  G,  corresponding  to 
the  two  hyperfine  lines.  Each  hyperfine  line  shows  a  single  positive  peak, 
indicating  that  the  magnetization  is  indeed  modulated  by  the  adiabatic  in¬ 
version.  A  small  peak  between  the  two  main  peaks  is  also  visible  and  is  due 
to  pairs  of  dopant  atoms  [29]. 
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Figure  3.7:  Force-detected  EPR  spectrum  of  the  low-concentration  sample 
for  interrupted  cyclic  adiabatic  inversion  for  T  =  6.5  K  and  CI/2it  =  10  MHz. 


Chapter  4 

E'  Centers  in  Si02 


This  chapter  describes  force-detected  EPR  experiments  investigating  singly- 
occupied  dangling  bond  defects  {E*  centers)  in  vitreous  silica.  This  type 
of  sample  is  important  because  similar  dangling  bond  defects  can  act  as 
trapping  sites  in  silicon  devices.  Furthermore  this  spin  system  is  the  leading 
candidate  for  our  initial  single  spin  detection  experiments  because  the  spin- 
lattice  relaxation  time  can  be  very  long  at  low  temperatures  {>  1  second), 
hopefully  allowing  cyclic  adiabatic  inversion  to  be  performed  coherently  for 
times  long  enough  for  signal  measurement. 

In  previous  EPR-MRFM  experiments [2,  15, 16,  30,  31,  32},  including  the 
one  described  in  the  previous  chapter,  modulation  of  the  sample  magneti¬ 
zation  has  been  accomplished  by  cyclic  saturation  using  modulation  of  the 
microwave  magnetic  field  Bi  and/or  modulation  of  the  polarizing  field  J5o[2].^ 
Cyclic  saturation  implicitly  depends  on  having  an  ensemble  of  spins  whose 
relaxation  time  is  much  shorter  than  the  modulation  period,  allowing  the 
magnetization  to  saturate  and  recover  nearly  instantaneously  as  the  sample 

^The  exception  is  the  use  of  the  interirupted  cyclic  adiabatic  inversion  described  in 
Section  3.2.4. 
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is  driven  in  and  out  of  resonance.  Since  the  concept  of  saturation  relies  on 
the  behavior  of  an  ensemble,  cyclic  saturation  is  not  a  suitable  technique 
for  future  MRFM  experiments  directed  towards  the  detection  of  individual 
electron  spins. 

Cyclic  adiabatic  inversion,  on  the  other  hand,  should  be  extendable  to 
the  single-spin  level.  In  cyclic  adiabatic  inversion,  the  spin  magnetization 
is  aligned  either  parallel  or  anti-parallel  to  the  effective  field  in  the  rotating 
frame  while  the  direction  of  the  effective  field  is  periodically  inverted  by 
modulating  the  frequency  of  the  Bi  field.  For  cyclic  adiabatic  inversion  to 
be  effective,  the  rotating  frame  relaxation  time[29]  must  be  long  enough  for 
a  large  number  of  inversion  cycles  to  be  completed,  allowing  a  measurable 
cantilever  vibration  to  build  up. 

Achieving  a  sufficiently  long  relaxation  time  is  problematic  for  electron 
spins  at  room  temperature  since  relaxation  times  are  commonly  in  the  mi¬ 
crosecond  or  nanosecond  range.  At  low  temperature,  however,  relaxation 
times  are  much  longer  and  a  number  of  electron  spin  systems  exhibit  suf¬ 
ficiently  slow  relaxation.  We  have  chosen  to  investigate  a  particularly  at¬ 
tractive  system:  E'  centers  in  vitreous  silica.  E'  centers  are  defect  sites 
characterized  by  singly  occupied  silicon  orbitals  [33].  For  spin  densities  on 
the  order  of  10^^  cm“®,  E'  centers  are  known  to  have  Ti  relaxation  times 
on  the  order  of  seconds  at  temperatures  as  high  as  5-10  K[34].  Such  long 
relaxation  times  are  well  suited  for  cyclic  adiabatic  inversion  experiments. 

The  work  in  this  chapter  serves  several  purposes.  First,  we  demonstrate 
that  cyclic  adiabatic  inversion  can  indeed  be  used  with  an  electron  spin  sys¬ 
tem  to  generate  a  detectable  oscillating  force.  After  describing  basic  MRFM 
signals  from  E'  centers,  we  show  measurements  of  non-equilibrium  magneti¬ 
zation  states  due  to  nutation  and  spin  echoes.  Finally,  we  study  the  relax¬ 
ation  behavior  of  E'  centers  during  cyclic  adiabatic  inversion  and  measure 
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the  rate  of  magnetization  decay  as  a  function  of  the  applied  Bi  field.  While 
the  experiments  presented  here  were  performed  using  a  micron-size  sample 
containing  millions  of  spins,  we  believe  that  the  magnetization  decay  results, 
in  particular,  will  be  useful  for  planning  future  single  spin  experiments.  Re¬ 
lated  experiments  on  bulk  samples  have  been  carried  out  using  a  conventional 
EPR  spectrometer  modified  to  detect  adiabatic  rapid  passage  signals[35,  36]. 

4.1  Force  detected  EPR 

4.1.1  Experiment 

The  experiments  were  performed  using  the  low-temperature  magnetic  reso¬ 
nance  force  detection  apparatus  used  for  the  experiments  described  in  the  pre¬ 
vious  chapter.  The  sample  material  was  synthetic  vitreous  silica  (Spectrosil) 
that  was  irradiated  by  ®°Co  7-rays  (4.4  x  10®  rad)  to  produce  E'  centers.^  The 
spin  concentration  was  measured  to  be  3.5  x  10^® /cm®  using  a  standard  spin¬ 
counting  technique[35,  36].  The  spin-lattice  relaxation  time  Ti  was  1.5  s  at 
5  K.  A  sample  piece  with  approximate  dimensions  50/xm  x50/.tm  x5/im  was 
glued  onto  a  commercial  single-crystal-silicon  cantilever[10]  with  a  spring  con¬ 
stant  k  ~  0.07  N/m  and  a  mechanical-resonant  frequency  fc  —  5.46  kHz  (with 
the  sample  loaded).  Angstrom-scale  motion  of  the  cantilever  was  monitored 
by  a  fiber-optic  interferometer[9]  and  synchronously  detected  with  a  lock-in 
amplifier.  In  order  to  speed  up  the  cantilever  response  time,  damping  feed¬ 
back  was  applied  to  the  cantilever  using  a  piezoelectric  element[ll,  12,  13]. 

The  sample  was  positioned  200  /im  above  the  center  of  a  microstripline 
resonator[19,  20]  that  produced  a  microwave  field  Bi  at  a;/27r  =  12.6  GHz. 
The  input  power  to  the  resonator  ranged  from  4  to  140  mW,  corresponding 


^Thanks  to  D.  L.  Griscom  for  generously  providing  this  material. 
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to  a  Bi  strength  of  0.23  to  1.35  G.  The  polarizing  field  Bq  was  a  combination 
of  the  homogeneous  field  from  a  superconducting  magnet  and  the  inhomo¬ 
geneous  field  (~  2300  G)  produced  by  a  permanent  magnet  located  800  fj,m 
away  from  the  sample.  The  permanent  magnet  created  a  field  gradient  of  2 
G//xm  and  served  to  generate  a  magnetic  force  on  the  sample.  All  experi¬ 
ments  were  performed  at  low  temperature  (5-10  K). 

4.1.2  Cyclic  adiabatic  inversion 

An  oscillatory  magnetic  force  from  the  resonant  spins  is  generated  by  a 
“probe  sequence” ,  which  consists  of  frequency  modulation  of  the  microwave 
field  such  that  a;(t)  =  wq  +  Aw(t),  where  coq  =  jBo-  Starting  with  Aw(0)  = 
— 27r  X  10  MHz  =  — the  microwave  field  is  switched  on  and  the  frequency  is 
modulated  such  that  A(o(t)  =  — J2cos{27r/ct).  After  several  thousand  cycles 
(t  =  0.4  s),  the  microwave  field  is  switched  off.  Provided  that  the  adia¬ 
batic  condition  (2.21)  is  satisfied  during  the  probe  sequence,  this  sequence 
induces  cyclic  inversion  of  the  resonant  magnetization  M,  which  follows  the 
direction  of  the  rotating  frcune  effective  field  as  described  in  Section  2.5.2. 
Consequently,  the  longitudinal  component  of  M  oscillates  at  a  frequency  fc, 
causing  the  cantilever  to  vibrate. 

Figure  4.1(a)  shows  the  cantilever- vibration  amplitude  during  a  typical 
probe  sequence  obtained  at  5.2  K  using  a  microwave  input  power  of  4  mW 
(Bi  =  0.23  G).  Damping  feedback  was  set  to  give  a  cantilever  response 
time  of  approximately  25  ms  and  a  lock-in  time  constant  of  3  ms  was  used. 
The  signal  was  averaged  over  roughly  30  probe  sequences,  with  4  seconds 
between  successive  sequences  to  allow  the  spins  to  repolarize.  As  shown  in 
Fig.  4.1(a),  the  cantilever  vibrated  in  response  to  the  probe  sequence  and 
a  peak  vibration  amplitude  of  ~  1  A-rms  was  achieved  approximately  25 


4.1.  FORCE  DETECTED  EPR 


39 


Figure  4.1:  (a)  Time  dependence  of  the  cantilever- vibration  amplitude.  The 
start  and  stop  of  the  probe  sequence  are  indicated  by  arrows.  The  probe 
sequence  consists  of  frequency  modulation  to  the  microwave  field  such  that 
Aa;(t)  =  —flcos(2irfct),  inducing  cyclic  adiabatic  inversion  of  the  electron 
spin  magnetization.  The  peak  frequency  deviation  Q/2'!r  was  10  MHz.  (b) 
Inverted  cantilever  response  obtained  when  sample  magnetization  was  adia- 
batically  inverted  by  a  frequency  ramp  before  the  probe  sequence. 
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ms  after  the  start  of  the  sequence.  After  the  initial  peak,  the  amplitude 
decayed  exponentially  with  a  decay  time  constant  of  approximately  55  ms. 
This  indicates  that  M  decayed  in  magnitude  while  it  was  modulated  by  the 
probe  sequence.  Taking  into  account  the  cantilever  response  time  and  lock-in 
time  constant,  we  estimate  the  decay  time  constant  of  the  magnetization  to 
be  25  ms.  This  magnetization  decay  will  be  discussed  further  in  Section  4.3. 

In  order  to  verify  that  frequency  modulation  of  Bi  indeed  causes  adiabatic 
inversion  of  M,  we  applied  a  single  additional  frequency  ramp  prior  to  the 
probe  sequence.  The  microwave  field  was  switched  on  with  Aco  =  Q  and 
then  ramped  to  Aa>  =  —Q.  This  ramp  should  adiabatically  invert  M.  The 
same  probe  sequence  as  before  was  then  applied:  Au;(t)  =  —Qcos(2irfct). 
As  shown  in  Fig.  4.1(b),  cantilever  vibration  was  excited,  but  with  opposite 
phase  compared  to  the  case  in  which  only  the  probe  sequence  was  applied. 
This  indicates  that  the  magnetization  before  the  start  of  the  probe  sequence 
was  indeed  adiabatically  inverted  by  the  preparatory  frequency  ramp. 

An  EPR  spectrum  can  be  obtained  by  measuring  the  cantilever  response 
to  the  probe  sequence  while  sweeping  the  field  from  the  superconducting 
magnet.  Figure  4.2  shows  the  peak  amplitude  of  the  cantilever  vibration  as 
a  function  of  Bq.  The  lineshape  is  symmetric  with  a  width  of  12  G.  This 
width,  considerably  broader  than  the  intrinsic  linewidth  for  E'  centers[37], 
results  from  the  use  of  the  large  gradient  with  a  sample  of  finite  thickness 
and  from  the  large  FM  deviation  f2/27r  =  20  MHz  (fi/7  —  7.1  G). 

4.1.3  Steady-state  magnetization  during  cyclic  adia¬ 
batic  inversion 

While  the  peak  amplitude  of  the  cantilever  vibration  reflects  the  state  of 
the  magnetization  just  prior  to  the  probe  sequence,  it  is  also  interesting 
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Figure  4.2:  Force-detected  EPR  spectrum.  The  peak  amplitude  of  the  can¬ 
tilever  vibration  is  plotted  as  a  function  of  the  polarizing  field  strength.  The 
peak  deviation  f2/27r  was  20  MHz,  the  cantilever  response  time  was  ~  10  ms, 
and  the  signal  was  averaged  over  50  scans.  The  peak  appears  when  Bq  =0)0/7 
at  the  sample  position. 
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to  investigate  the  steady-state  cantilever  vibration  achieved  at  the  end  of 
the  probe  sequence.  In  Fig.  4.1,  the  cantilever-vibration  amplitude  decays 
almost  to  zero  by  the  end  of  the  probe  sequence;  however,  we  found  that  this 
is  not  always  the  case  and  that  the  steady-state  amplitude  depends  on  the 
strengths  of  both  Bq  and  Bi.  Figure  4.3(a)  shows  the  cantilever- vibration 
amplitude  during  the  probe  sequence  with  a  microwave  input  power  of  52  mW 
{Bi  =  0.82  G).  After  the  initial  excitation,  the  vibration  amplitude  decays 
and  reaches  a  steady-state  value  of  approximately  0.3  A-rms. 

The  steady-state  value  of  cantilever-vibration  amplitude  changes  sign 
as  the  external  field  is  swept,  leading  to  the  bipolar  spectrum^  shown  in 
Fig.  4.3(b).  The  cantilever- vibration  amplitude  reaches  a  positive  steady- 
state  value  when  Bq  is  slightly  above  resonance  {Bq  —  cj/'y  >  0),  and  a 
negative  value  when  Bq  is  below  resonance  [Bq  —  u/j  <  0).  These  results 
can  be  understood  in  terms  of  spin  dynamics  in  the  rotating  frame.  In  the 
steady  state,  M  is  locked  parallel  to  Beir  if  Bq  —  0^/7  >  0  and  antipaxallel  to 
Beff  if  Bo  —  a; /7  <  0[29].  When  cyclically  modulated  by  the  FM  sequence, 
the  magnetization  antiparallel  to  Beff  generates  an  oscillatory  force  that  is 
180°  out  of  phase  with  respect  to  that  generated  when  M  is  parallel  to  Beff. 

4.2  Nonequilibrium  magnetization  measure¬ 
ments 

4.2.1  Nutation 

The  initial  peak  amplitude  of  the  cantilever  vibration  reflects  the  state  of 
longitudinal  magnetization  at  the  start  of  the  probe  sequence.  This  fact  can 

®The  smeill  background  signal  of  ~  0.02  A  observed  for  field  values  far  from  resonance 
is  due  to  spurious  microwave-feedthrough  effects. 
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(a) 


Figure  4.3:  (a)  Time  dependence  of  the  cantilever-vibration  amplitude  with 
a  microwave  input  power  of  52  mW  {Bi  =  0.82  G).  The  vibration  amplitude 
reaches  a  steady-state  value  by  the  end  of  the  probe  sequence,  (b)  The 
steady-state  cantilever-vibration  amplitude  as  a  function  of  the  Bq  strength. 
Absolute  vibration  amplitude  is  smaller  than  in  (a)  because  greater  damping 
feedback  was  used. 
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be  used  to  probe  nonequilibrium  states  of  magnetization[28].  Here,  we  show 
nutation  and  spin  echo  experiments. 

Nutations  can  be  measured  by  applying  a  microwave  pulse  of  width  tp 
before  the  probe  sequence.  The  pulse  causes  resonant  spins  to  precess  about 
Beff  =  in  the  rotating  frame  with  an  angular  frequency  7H1.  As  a  result, 
the  longitudinal  component  of  M  is  given  by  Mz  =  Mq  cos['yBitp),  where  Mq 
is  the  thermal  equilibrium  magnetization.  The  probe  sequence  following  the 
pulse  generates  an  oscillatory  force  proportional  to  Mz  that  excites  the  can¬ 
tilever  vibration.  Figure  4.4  shows  the  peak  cantilever-vibration  amplitude 
as  a  function  of  the  microwave  pulse  width  tp.  A  microwave  input  power  of 
130  mW  was  used  for  the  pulse  and  4  mW  was  used  for  the  probe  sequence. 
The  signals  show  an  oscillation  with  a  270  ns  period  and  gradually  decreasing 
amplitude.  The  oscillation  is  superimposed  on  a  positive  offset.  Since  the 
oscillation  period  corresponds  to  2‘Kf{'yBi),  the  amplitude  of  Bi  for  the  130 
mW  microwave  pulse  is  determined  to  be  1.3  G.  This  result  calibrates  the  Bi 
field  at  the  location  of  the  sample  generated  by  the  microstripline  resonator. 

The  positive  offset  is  due  to  off-resonance  spins  in  the  sample  (i.c.,  due 
to  the  field  gradient  and  the  finite  sample  thickness).  With  the  frequency 
deviation  n/27r  =  10  MHz  used  for  the  probe  sequence,  spins  that  are 
off-resonance  by  |AB|  =  J2/7  ~  3.6  G  still  contribute  to  the  force  signal 
significantly [17,  28].^  Rather  than  precessing  about  J5iX,  these  off-resonance 
spins  will  precess  in  a  cone  about  Beff  =  HiX-1- AHz,  which  is  tilted  70°  from 
the  X  axis  (for  Bi  =  1.3  G).  Spins  precessing  in  this  cone  will  always  have  a 
positive  z  component,  leading  to  a  positive  offset  in  the  nutation. 

The  decay  of  the  nutation  is  also  attributed  to  off-resonance  spins.  For  the 


^See  Section  ??. 
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Figure  4.4;  Force-detected  nutation  signal.  The  peak  amplitude  of  the  can¬ 
tilever  vibration  is  plotted  as  a  function  of  the  width  of  a  microwave  pulse 
tp  applied  prior  to  the  probe  sequence.  A  microwave  power  of  130  mW  was 
used  for  the  pulse,  and  4  mW  was  used  for  the  probe  sequence.  The  dotted 
line  is  a  fit  to  a  damped  cosinusoidal  oscillation  with  an  offset. 
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ofF-resonance  case,  the  precession  frequency  7|Beff|  differs  from  7B1,  causing 
spins  in  different  regions  of  the  sample  to  process  at  different  rates.  The 
interference  caused  by  a  range  of  precession  frequencies  leads  to  a  gradual 
decay  of  the  oscillation  amplitude  in  the  total  nutation  signal.  This  effect 
dominates  over  the  decay  due  to  spin-spin  relaxation  processes  {T2  decay), 
which  are  related  to  the  much  smaller  homogeneous  broadening. 


4.2.2  Spin  echoes 

Spin  echoes  play  an  important  role  in  many  magnetic  resonance  techniques. 

In  conventional  magnetic  resonance,  a  7r/2-r-7r  pulse  sequence  is  used  to 
generate  an  echo.  After  the  7r/2  pulse,  the  spread  in  precession  rates  of 
the  magnetization  due  to  local  fields  and  the  inhomogeneity  of  Bq  causes  a 
reversible  decay  of  transverse  magnetization.  When  a  it  pulse  is  applied  after 
a  time  r,  spins  are  refocused  at  t  =  2t,  leading  to  an  echo,  the  size  of  which 
is  given  by  Moexp{—2TjT2).  Since  the  MRFM  technique  described  here  is 
sensitive  to  the  longitudinal  component  of  M,  the  spin  echo  is  detected  by 
applying  a  7r/2-r-7r-T'-7r/2  pulse  sequence.  The  last  it j 2  pulse  rotates  the 
refocused  magnetization  to  the  longitudinal  direction,  which  is  mechanically 
detected  by  a  subsequent  probe  sequence. 

The  experiment  was  carried  out  with  the  pulse  sequence  shown  in  Fig.  4.5(a). 
The  peak  amplitude  of  the  cantilever  vibration  was  measured  while  changing 
the  time  interval  between  the  first  7r/2  pulse  and  the  tt  pulse  and  keeping 
a  fixed  time  interval  tb  between  the  two  7r/2  pulses.  This  has  the  effect  of 
“walking”  the  echo  through  the  last  it  j 2  pulse.  The  widths  of  the  7r/2  and 
TT  pulses  were  70  ns  and  140  ns,  respectively.  Figure  4.5(b)  shows  the  result 
for  Tb  =  1  /Its.  A  peak  is  observed  at  ta  —  500  ns  =  tb/2,  corresponding  to 
the  refocus  condition  for  a  spin  echo. 
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Figure  4.5:  (a)  Pulse  sequence  used  for  the  spin  echo  experiment.  A  mi¬ 
crowave  power  of  130  mW  was  used  for  the  pulses,  and  4  mW  was  used  for 
the  probe  sequence.  Widths  of  7r/2  and  tt  pulses  were  70  ns  and  140  ns,  re¬ 
spectively.  (b)  The  peak  amplitude  of  the  cantilever  vibration  as  a  function 
of  ta,  with  a  fixed  tb  =  1  /is.  A  peak  appears  when  ta  —  Tb/2,  correspond¬ 
ing  to  a  spin  echo,  (c)  Spin  echoes  observed  with  various  tb.  The  amplitude 
of  the  spin  echoes  gradually  decreases  as  tb  increases. 
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The  results  for  various  tb  are  shown  in  Fig.  4.5(c).  In  each  case,  a  peak 
appears  when  ta  —  representing  spin  echoes  at  various  time  intervals. 

With  increasing  tb,  the  peak  height  gradually  decreases,  presumably  due  to 
the  T2  decay  as  in  conventional  spin  echo  experiments  [38].  The  observed 
decay  time  indicates  a  T2  on  the  order  of  5  fis.  For  comparison,  T2  of  2.1 
ixs  has  been  measured  previously  using  conventional  techniques  on  a  silica 
sample  containing  a  slightly  higher  concentration  of  spins[39]. 
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4.3  Magnetization  decay  during  adiabatic  in¬ 
version  cycles 

One  key  issue  for  the  force  detection  of  electron  spins  is  the  relaxation  rate 
of  the  spins  while  they  are  subjected  to  cyclic  adiabatic  inversion.  As  men¬ 
tioned  earlier,  the  magnetization  decays  during  the  force  detection  sequence 
(Fig.  4.1).  For  a  successful  MRFM  experiment,  the  spin  magnetization  must 
last  long  enough  to  excite  the  cantilever  vibration  to  a  detectable  level.  The 
decay  of  magnetization  can  be  due  to  a  variety  of  effects,  including  spin- 
lattice  relaxation  in  the  rotating  frame,  violation  of  the  adiabatic  condition, 
and  spin-spin  relaxation  processes. 

As  shown  in  Fig.  4.6,  the  spin-lattice  relaxation  time  in  the  rotating 
frame  Tip  can  be  measured  by  staying  spin-locked  on  resonance  for  a  variable 
time  Tiock  before  the  probe  sequence.  The  data  in  Fig.  4.6  shows  the  peak 
cantilever- vibration  amplitude  as  a  function  of  the  spin-locking  time  for  Bi  = 
0.23  G.  The  signal  decreases  exponentially  towards  a  steady-state  value  as 
Tlock  increases.  The  decay  time  constant  was  estimated  to  be  Tip  ~  190 
ms,  in  reasonable  agreement  with  the  value  measured  using  conventional 
techniques  [35,  36]. 

The  observed  magnetization  decay  time  constant  during  repeated  adia¬ 
batic  reversals  was  found  to  be  significantly  shorter  than  Tip.  For  example, 
the  observed  decay  time  constant  in  Fig.  4.1  was  only  25  ms.  To  further  in¬ 
vestigate  this  effect,  we  used  the  trapezoidal  FM  waveform  shown  in  Fig.  4.7, 
instead  of  the  usual  sinusoidal  frequency  modulation.  The  slope  of  the  lineaj 
frequency  ramp  in  the  waveform  was  set  such  that  \'y~^du>/dt\  =  \dBes/dt\  = 
3.0  X  10®  G/s. 

The  frequency  modulation  of  Bi  causes  an  oscillation  of  longitudinal  mag- 
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Figure  4.6:  The  peak  amplitude  of  the  cantilever  vibration  as  a  function  of 
the  spin-lock  time  riock  before  the  probe  sequence.  The  dotted  line  is  a  fit  to 
an  exponential  decay  with  an  offset.  The  signal  decays  with  a  time  constant 
Tip  Ci  190  ms.  The  inset  shows  the  FM  waveform  used  for  the  measurement. 
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Figure  4.7:  Trapezoidal  FM  waveform  used  to  study  magnetization  decay. 
The  slope  of  the  trapezoid  is  set  such  that  \'y~^duj/dt\  =  |dBeff/dt|  =  3.0x  10® 
G/s. 


52 


CHAPTER  4.  E'  CENTERS  IN  SIO2 


netization  such  that 


M,{t)  =  M{t) 


Aoj(t) 


(4.1) 


where  M (t)  is  the  magnitude  of  the  magnetization,  which  varies  slowly  com¬ 
pared  to  the  5.5  kHz  modulation  frequency.  The  first-harmonic  Fourier 
component  (in  rms)  of  M2(t)  is  approximately  given  by  {2\/2jTr)M{t)  for 
a  Au}{t)  with  a  trapezoidal  waveform.  If  M{t)  slowly  decays  such  that 
M{t)  =  Moexp{—tfTjn),  the  rms  amplitude  of  the  resulting  oscillatory  force 
is  given  by 

2V^ 


F{t)  =  ^—^MoGexp  ( — +  Fb, 

TT  \  Tm  / 


(4.2) 


where  G  is  the  field  gradient  and  Fb  is  a  small  background  oscillatory  force 
that  is  included  to  account  for  spurious  feedthrough  of  the  microwave  mod¬ 
ulation. 

The  rms  amplitude  of  the  cantilever  response  A{t)  is  given  by  the  convo¬ 
lution 

A{t)  =  r  h{t')F{t  -  t')dt\  (4.3) 

Jo 

where 

(-0  (-) 


h{t)  =  exp 


kTc 

is  the  impulse  response  function  of  the  cantilever.  Here,  Qeff  and  Tg  = 
Qeff/7r/c  are  the  effective  Q  factor  and  1/e  response  time  of  the  cantilever, 
respectively,  taking  into  account  the  damping  feedback.  Equations  (4.2)  and 
(4.3)  lead  to  the  cantilever  response 


m 


QeS  f  2\/2  MqG 
A:  \  TT  l-rc/r„ 


K  (■;■)  “  (“0] + 

(4.5) 

A  typical  cantilever  response  for  the  trapezoidal  modulation  sequence  is 
shown  in  Fig.  4.8.  The  observed  cantilever  response  was  fit  to  Eq.  (4.5)  using 
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Figure  4.8:  A  typical  cantilever  response  for  the  trapezoidal  modulation.  The 
line  is  a  fit  using  the  function  in  Eq.  (4.5).  The  magnetization  decay  time 
Tm  is  found  to  be  46  ms  in  this  case.  Bi  —  1.35  G. 
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three  adjustable  parameters  Mq,  Tm,  and  Fb.  The  other  parameters  were 
determined  independently.  For  Fig.  4.8,  we  found  that  Mq  =  2.4  x  10~^^ 
emu  =  2.6  x  where  (ib  is  the  Bohr  magneton,  and  Tm  =  46  ms.  The 

uncertainty  in  these  results  is  estimated  to  be  on  the  order  of  10%. 

Magnetization  decay  time  constants  Tm  obtained  for  various  Bi  strengths 
are  plotted  in  Fig.  4.9.  Values  of  adiabaticity  parameter  'yBl/\dB^([/dt\, 
defined  as  the  ratio  of  the  right  hand  side  and  left  hand  side  of  (2.21), 
axe  also  indicated  in  the  figure.  It  was  found  that  Tm  increases  as  Bi  in¬ 
creases,  suggesting  that  the  magnetization  decay  is  due  to  the  violation  of 
the  adiabatic  condition,  \dB^{{/dt\  <  'fB^.  At  the  maximum  Bi,  where 
jBl/\dBes/dt\  ~  10,  Tm  =  46  ms  and  it  appears  to  be  reaching  a  plateau, 
indicating  that  loss  due  to  non-adiabaticity  is  no  longer  dominant.  Never¬ 
theless,  Tm  is  still  significantly  shorter  than  Tip  ~  190  ms.  An  additional 
mechanism  for  magnetization  decay,  such  as  rotating  frame  cross-relaxation 
processes  among  spins  in  inhomogeneously  broadened  lines,  might  be  respon¬ 
sible  for  the  relatively  short  36]. 

The  magnetization  decay  time  during  cyclic  adiabatic  inversion  may  have 
implications  for  future  single  spin  experiments  since  the  spin  lifetime  deter¬ 
mines  the  detection  bandwidth  that  must  be  used.  The  spin  lifetime  of  ~  50 
ms  implies  the  detection  bandwidth  of  ~  20  Hz.  A  sharp  magnetic  tip  (50 
um  radius)  generates  a  field  gradient  on  the  order  of  10  G/A  at  a  position 
~  150  A  below  the  tip.  For  an  electron  spin  located  at  this  position  and  un¬ 
dergoing  cyclic  adiabatic  inversion,  the  resulting  interaction  force  is  7  X  10 
N-rms.  To  detect  such  a  small  force  in  a  detection  bandwidth  of  20  Hz,  a 
force  resolution  of  2  x  10“^^  N/\/Hz  or  better  is  needed  to  achieve  unity 
signal-to-noise  ratio.  Such  force  resolution  was  recently  demonstrated  using 
ultrathin  single-crystal  silicon  cantilevers[40,  41].  Assuming  that  the  spin  re¬ 
laxation  time  is  not  adversely  affected  by  the  close  proximity  of  the  magnetic 


Magnetization  decay  time  (ms) 


y  B,"/ldB  Jdtl  =  1 
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Figure  4.9:  Magnetization  decay  time  as  a  function  of  Bi.  Arrows  indicate 
values  of  adiabaticity  parameter  'yBi/\d'Be{i/dt\. 
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tip,  the  long  spin  lifetime  of  E'  centers  make  them  a  suitable  candidate  for 
future  single  spin  MRFM  experiments. 


Chapter  5 

EPR  and  FMR  imaging 


Magnetic  resonance  force  microscopy  (MRFM)  has  previously  been  per¬ 
formed  with  the  sample  mounted  on  a  micromechanical  cantilever  positioned 
close  to  a  millimeter-size  field  gradient  source.  In  this  chapter,  we  describe 
an  improved  microscope  configuration,  where  the  gradient  source  consists  of 
a  magnetic  tip  mounted  directly  on  the  cantilever.  This  allows  a  much  wider 
variety  of  samples  to  be  investigated  and  greatly  improves  the  convenience 
of  the  technique.  We  have  used  this  improved  microscope  to  compare  and 
contrast  two-  and  three-dimensional  imaging  characteristics  of  paramagnetic 
and  ferromagnetic  samples. 


5.1  Improved  MRFM  setup 

The  design  of  the  microscope  is  shown  schematically  in  Figs.  5.1,  5.2,  and  5.3. 
The  entire  system  is  suspended  in  a  vacuum  can  with  Teflon  support  frames 
to  ensure  adequate  vibration  isolation.  The  vacuum  can  fits  into  a  liquid- 
helium  dewer  with  a  superconducting  magnet.  The  microscope  consists  of 
a  gradient  source-cantilever  combination  and  a  sample  that  can  be  scanned 
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in  three  dimensions.  The  gradient  source  is  a  magnetic  tip  (typically  a  10- 
100  micron-size  sphere)  glued  onto  the  cantilever,  the  motion  of  which  is 
measured  by  a  fiber-optic  interferometer [9].  The  sample  is  placed  on  top 
of  a  microstripline  resonator [19,  20]  that  generates  the  12-GHz  microwave 
field  Bi  used  to  induce  magnetic  resonance.  The  sample  and  resonator  are 
mounted  on  a  piezoelectric  tube  scanner. 

In  this  design,  the  spatial  resolution  of  MRFM  is  aimed  at  a  micron  level, 
requiring  a  total  scan  range  on  the  order  of  100  fiva  at  low  temperature. 
We  used  a  15-cm-long  piezoelectric  tube  scanner  to  achieve  this.  Two  tube 
scanners  (3  inch  long  x  0.25  inch  diameter  x  0.02  inch  wall-thickness) [42] 
were  joined  together  to  achieve  a  scan  range  of  approximately  100  fxm  lat¬ 
erally  (x  and  y)  and  5  //m  in  the  z  direction  at  low  temperature  with  a 
maximum  applied  voltage  of  il50  V.^  The  piezoelectric  tube  is  surrounded 
by  an  outer  support  tube  attached  to  a  tripod  assembly,  which  is  held  with 
three  springs  against  three  adjustment  screws.  One  of  these  screws  is  coupled 
to  an  external  shaft  for  coaxse  adjustment  of  the  tip-sample  distance  during 
operation[8]. 

In  order  to  provide  thermal  conduction  from  liquid  helium,  most  of  the 
instrument,  including  the  sample/resonator  holder,  the  outer  support  tube 
of  the  scanner,  and  the  side  frame  of  the  microscope,  is  made  of  oxygen-free 
high-conductivity  (OFHC)  copper.  The  copper  braids,  attached  by  soldering, 
connect  to  the  bottom  flange  of  the  vacuum  can  immersed  in  liquid  helium. 
Since  the  copper  braid  and  a  semi-rigid  coax  to  provide  the  microwave  to  the 
resonator  are  attached  to  the  sample/resonator  holder  on  top  of  the  15-cm- 
long  tube  scanner,  special  attention  is  required  in  handling  the  microscope 
not  to  break  the  rather  fragile  tube  scanner.  One  disadvantage  of  using  such 

^The  scan  range  was  calibrated  by  imaging  a  calibration  grid  using  contact-mode  atomic 
force  microscopy. 
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Figure  5.1:  Schematic  diagram  of  the  improved  low-temperature  magnetic 
resonance  force  microscope.  The  microscope  is  suspended  in  a  vacuum  can 
that  fits  into  a  liquid-helium  dewer  with  a  superconducting  magnet. 
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a  long  tube  scanner  is  that  the  positioning  of  the  sample  holder  is  susceptible 
to  environmental  vibration  coupled  either  directly  via  the  support  structure 
or  through  the  coax  and  copper  braid.  The  stability  of  the  sample  holder 
was  estimated  to  be  on  the  order  of  0. 1-1.0  /im  in  x  and  y,  and  0.1  /um  in  z. 

The  cantilever  chip  is  glued  onto  the  end  of  a  stiff  stainless  steel  flexure 
beam  that  is  deflected  by  a  fine-pitch  adjustment  screw  as  shown  in  Fig.  5.2. 
The  screw  is  coupled  to  a  second  external  shaft.  Between  the  screw  and 
the  flexure  beam  is  a  piezoelectric  disk  used  to  fine-adjust  the  z  position 
of  the  cantilever  relative  to  the  fiber  end.  The  piezoelectric  disk  is  also 
used  to  drive  the  vibration  of  the  cantilever  and  provide  damping  feedback 
in  some  casesfll,  12,  13].  The  x-y  position  of  the  fiber  end  is  adjusted  at 
room  temperature  with  the  positioner  consisting  of  a  cylindrical  flexure  and 
adjustment  screws. 

Magnetic  resonance  in  the  sample  is  induced  by  the  microwave  field  Bi 
from  the  resonator  in  the  presence  of  the  static  field  Bq,  which  is  a  combina¬ 
tion  of  the  homogeneous  field  B^xt  applied  in  the  z  direction  by  an  external 
superconducting  magnet  and  the  inhomogeneous  field  Btip  from  the  mag¬ 
netic  tip.  Bcxt  also  served  to  magnetize  the  magnetic  tip  to  its  saturation 
magnetization. 

To  generate  a  kilohertz-frequency  oscillatory  magnetic  force  from  the  on- 
resonance  spins,  Bi  was  frequency  modulated  so  as  to  cyclically  saturate  the 
spin  magnetization.  The  modulation  frequency  cOm  was  chosen  to  be  equal 
to  the  mechanical  resonant  frequency  of  the  cantilever  tOc.  This  results  in  a 
bipolar  spectral  lineshape  of  force-detected  magnetic  resonance  as  explained 
in  Chapter  3.  The  force-detected  magnetic  resonance  signal,  measured  as  a 
vibration  amplitude  of  the  cantilever,  is  recorded  while  the  sample  is  scanned 
in  x-y  with  various  values  of  B^xt- 

The  presence  of  the  tip  on  the  cantilever  results  in  some  unwanted  side 
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Figure  5.2:  Close-up  view  of  the  microscope  from  the  side. 
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Figure  5.3:  Schematic  setup  of  the  experiment 
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effects  observed  during  tbe  x-y-Bcxt  scan.  The  mechanical  resonant  fre¬ 
quency  of  the  cantilever  is  affected  by  tip-sample  interaction  forces  and  by 
the  force/torque  exerted  by  Therefore,  at  each  point  of  the  x-y-Bext 

scan,  u>c  is  first  measured  in  order  to  set  ojrn  before  detecting  the  magnetic 
resonance  force  signal.  To  ensure  proper  operation  of  the  fiber  interferom¬ 
eter,  the  cantilever-fiber  spacing  is  corrected  by  a  piezoelectric  actuator  at 
the  base  of  the  cantilever. 


Figure  5.4  shows  a  block  diagram  of  the  experiment.  A  personal  com¬ 
puter  (PC)  controls  the  sample  position  and  B^xt-  The  cantilever-fiber  spac¬ 
ing  is  controlled  by  servo  electronics  to  maintain  the  proper  operation  of  the 
interferometer.  At  each  point  of  the  scan,  a  calibration  sequence  for  the  can¬ 
tilever  frequency  is  first  initiated.  A  pulse  generator  (HP8118A)  sends  pulses 
to  open  the  gate  to  drive  the  cantilever  vibration  with  the  piezoelectric  disk, 
and  to  activate  a  phase-locked  loop  (PLL)  consisting  of  a  lock-in  amplifier 
(PAR124)  to  detect  the  phase  error,  which  is  fedback  to  the  frequency  control 
of  the  function  generator  (FG5010)  driving  the  cantilever  vibration.  After 
the  calibration  sequence,  the  frequency  modulation  of  the  microwave  genera¬ 
tor  (HP83732A)  is  turned  on  at  the  calibrated  cantilever  resonant  frequency. 
The  cantilever  vibration  amplitude  resulting  from  the  oscillatory  force  due  to 
magnetic  resonance  is  then  measured  by  another  lock-in  amplifier  (SR530) 
and  recorded  by  the  PC.  K  necessary,  damping  feedback  to  reduce  the  can¬ 
tilever  response  time  is  used  by  applying  a  90°-phase-shifted  signal  to  the 
piezoelectric  disk. 
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Figure  5.4:  Block  diagram  of  the  experiment. 
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5.2  Cantilever  with  a  magnetic  tip  in  the  ex¬ 
ternal  field 

As  a  magnetic  tip,  we  used  a  NdFeB  sphere  of  ~10  fim  radius,^  a  ferrite  par¬ 
ticle  (irregular  shape),  and  nickel  spheres  of  ~10  fim  and  ~55  fim  radii.^  The 
magnetic  tip  was  glued  onto  a  commercial  single-crystal-silicon  cantilever  [10]. 
Figure  5.5  shows  optical  micrographs  of  10  fim  and  55  ^m  radius  nickel 
spheres  mounted  on  the  cantilevers. 

The  behavior  of  the  cantilever  with  the  magnetic  tip  is  studied  by  measur¬ 
ing  the  dc  deflection  and  resonant  frequency  of  the  cantilever  in  the  applied 
external  field.  Figure  5.6  shows  the  dc  deflection  of  the  cantilever  with  a  10- 
;um-radius  nickel  sphere  (Fig.  5.5(a)).  The  deflection  increases  as  the  field  is 
increased  in  both  polarities,  with  greatest  effect  at  lower  field  (below  ~3000 
G). 

The  magnetic  force  acting  on  the  magnetic  tip  is 

P  =  MV—,  (5.1) 

where  M  is  the  magnetization  of  the  tip  assumed  to  be  parallel  to  the  field 
B  applied  in  the  z  direction  and  V  is  the  volume  of  the  tip.  The  cantilever  is 
not  positioned  exactly  at  the  center  of  the  superconducting  magnet,  resulting 
in  a  small  field  gradient  on  the  order  of  10  G/cm  for  B  —  1000  G.  For  a 
soft  ferromagnetic  sphere,  M  varies  linearly  with  B  up  to  where 

Ms  is  the  saturation  magnetization.  For  nickel,  Ms  ^  510  emu/cm®  at  low 
temperature.  Since  dBjdz  is  proportional  to  B,  the  linear  dependence  of  M 
on  B  up  to  AkMsIZ  ~  2140  G  leads  to  F  oc  B^.  In  higher  fields  where  M 


^Available  from  Particle  Technology,  Inc.,  Hanover,  Maryland. 
^Available  from  Alfa  Aester,  1-800-343-0660. 
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Figure  5.5:  (a)  Optical  micrographs  of  the  cantilever  with  a  nickel  sphere 
with  r  ~  10  /zm.  A  pyramidal-shape  silicon  tip  is  also  visible,  (b)  The 
cantilever  with  a  nickel  sphere  with  r  ~  55  fim. 
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reaches  a  constant  value  of  Ms,  F  is  expected  to  be  linear  with  B.  The  curve 
in  Fig.  5.6  agrees  qualitatively  with  the  expected  behavior. 

The  curve  shows  a  little  hysteresis,  indicating  the  presence  of  anisotropy 
presumably  due  to  imperfect  spherical  shape  or  crystalline  anisotropy.  If 
there  is  anisotropy,  a  torque  also  acts  on  the  tip  since  the  magnetization 
tends  to  align  with  the  easy  axis  as  well  as  the  external  field.  Let  6  be  the 
angle  between  the  easy  axis  and  the  z  direction  when  the  cantilever  is  not 
deflected,  a  be  the  angle  between  the  easy  axis  and  the  magnetization,  and 
ip  be  the  angle  of  the  cantilever  deflection  as  shown  in  Fig.  5.7.  In  the  case 
of  the  crystalline  anisotropy,  the  total  energy  can  be  modeled  as 


E  = 


dB 

K  sin^  a  —  BM  cos(^  —  a  —  (p)  + 


V  +  sin^  ip,  (5.2) 


where  K  is  the  anisotropy  factor,  h  and  I  axe  the  spring  constant  and  the 
length  of  the  cantilever,  respectively.  The  first  term  in  E  is  the  anisotropy 
energy,  the  second  and  the  third  terms  axe  the  energies  in  the  applied  field, 
and  the  last  term  is  the  potential  energy  of  the  cantilever.  The  cantilever 
deflection  can  be  determined  by  minimizing  E  with  respect  to  a  and  ip  with 
given  B  and  M,  which  is  assumed  to  be  proportional  to  B  up  to  AirMsjZ. 
Even  if  B  is  perfectly  uniform  {dBjdz  —  0),  this  leads  to  a  finite  deflection. 
For  small  B,  we  obtain  a  ~  0  and  ip  a  B"^.  For  large  B,  where  M  is  saturated, 
ip  reaches  a  constant  value  6.  The  dc  deflection  shown  in  Fig.  5.6  is  likely 
to  be  resulting  from  both  the  force  due  to  dBjdz  and  the  torque  due  to  the 
anisotropy. 

Figure  5.8  shows  the  change  of  the  resonant  frequency  as  a  function  of 
the  field  for  the  same  cantilever.  The  resonant  frequency  at  zero  field  is 
3.97  kHz.  The  frequency  increases  up  to  ~2100  G  and  then  decreases,  reach¬ 
ing  a  plateau  above  ~5000  G.  There  is  also  a  little  hysteresis  due  to  an 
apparent  ~600  G  of  coercivity. 
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Figure  5.6:  The  dc  deflection  of  the  cantilever  with  a  nickel  sphere  (r  ~ 
/Lim)  as  the  field  is  swept  forward  and  backward. 
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Figure  5.7:  Model  of  the  cantilever  deflection. 


The  magnetic  interaction  modifies  the  effective  spring  constant  of  the 
cantilever  according  to 

dF 

A.=  ^,  (6.3) 

which,  in  the  absence  of  the  anisotropy,  is  given  by 


^  =  mv^21. 
dz  dz^ 


(5,4) 


Near  the  center  of  the  superconducting  magnet,  which  can  be  assumed  to 
be  a  simple  solenoid,  this  term  has  the  effect  of  increasing  the  resonant 
frequency.  Since  both  M  and  d^Bldz^  are  proportional  to  B  for  low  field,  the 
frequency  change  is  expected  to  be  proportional  to  B^  up  to  Air M,  1 3  ~  2140 
G,  beyond  which  it  is  expected  to  be  proportional  to  B.  The  torque  due  to 
the  anisotropy  also  has  the  effect  of  increasing  the  frequency,  which  varies 
as  oc  B^  for  low  field  and  reaches  a  plateau  for  high  field.  The  result  shown 
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in  Fig.  5.8  agrees  with  the  expected  behavior  for  low  field.  The  frequency 
decrease  in  high  field,  however,  does  not  agree  with  the  expectation  based 
on  the  model  described  in  Eq.  (5.2),  indicating  an  additional  mechanism 
affecting  the  resonant  frequency  of  the  cantilever. 

The  presence  of  the  tip  on  the  cantilever  also  has  a  significant  effect  on 
the  quality  factor  of  the  cantilever  in  the  field.  Figure  5.9  shows  Q  for  the 
cantilevers  with  10  /im  and  55  fim  radius  nickel  spheres  (Fig.  5.5).  Q  is 
reduced  from  more  than  100  000  at  zero  field  as  the  field  is  increased.  The 
effect  is  more  significant  for  the  bigger  sphere,  for  which  Q  is  reduced  to 
below  100  at  B  =10  000  G.  This  effect  is  due  to  micromagnetic  loss  caused 
by  motion  of  the  tip  magnetization  when  the  cantilever  vibrates  in  the  field. 


5.3  EPR  of  phosphorus-doped  silicon 

As  described  in  Chapter  3,  electron  paxamagnetic  resonance  (EPR)  of  phosphorus- 
doped  silicon  was  studied  using  the  low-temperature  magnetic  resonance 
force  detection  apparatus.  It  was  found  that  force  detection  can  be  used  for 
spectroscopic  measurements  in  micron-size  samples,  and  the  hyperfine  split¬ 
ting  of  the  EPR  line  was  observed  for  the  low  dopant-concentration  sample 
(iVu  ~  8  X  10^®  cm-3). 

Here,  we  compare  imaging  characteristics  of  low-  and  high-concentration 
phosphorus-doped  silicon  using  the  improved  low-temperature  MRFM  appa¬ 
ratus. 

Samples  were  approximately  100-//m-thick  phosphorus-doped  silicon  pieces 
(4  mm  X  2  mm).  Two  samples  with  dopant  concentrations  Nd  ~  8  x  10^® 
cm“®  and  4  x  10^®  cm“®  were  used.  A  55-/im-radius  nickel  sphere  was  used 
as  a  magnetic  tip  (Fig.  5.5(b)). 
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Figure  5.8:  The  resonant  frequency  of  the  cantilever  with  a  nickel  sphere 
(r  ~  10  fj,m)  as  the  field  is  swept  forward  and  backward. 
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(a) 


B(G) 


Figure  5.9:  (a)  Q  of  the  cantilever  with  a  nickel  sphere  (r  ~  10  //m)  as  a 
function  of  the  field,  (b)  Q  of  the  cantilever  with  a  nickel  sphere  (r  ~  55 
fj,m)  as  a  function  of  the  field. 
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Since  the  samples  are  uniformly  doped  and  the  lateral  dimensions  are 
much  larger  than  the  tip  size,  scanning  in  x  and  y  does  not  provide  any 
information  on  the  imaging  characteristics.  Instead,  EPR  spectra  were  ob¬ 
tained  by  sweeping  the  external  field  B^xt  with  Afarious  tip-sample  distances. 


5.3.1  Low-dopant-concentration  silicon 

Figure  5.10  shows  EPR  spectra  for  low-concentration  sample  {Nd  ~  8  x 
10^®  cm“^).  Tip-sample  distances  were  approximately  80  /im  for  (a),  60  fxm 
for  (b),  and  40  /j.m  for  (c),  respectively.  A  bipolar  signal  was  observed  at 
~4300  G  independent  of  the  tip-sample  distance.  The  peah-to-peak  force 
was  approximately  100  fN. 

A  smaller  peak,  indicated  by  arrows,  was  also  observed  downfield  from 
the  bipolar  peak.  The  position  of  the  small  peak  moved  to  lower  field  as  the 
tip-sample  distance  was  decreased. 

Magnetic  resonance  is  induced  for  spins  within  the  resonant  slice,  where 
|Bext  +  Btipl  4300  G  for  electron  spins  in  the  12-GHz  microwave  field.  If 
|Bext|  ^  I  Btipl,  which  is  typically  the  case  in  the  experiments,  the  resonant 
slice  is  a  paraboloidal-shape  contour  of  the  z  component  of  Btip,  which  we 
designate  Figure  5.11  shows  contours  of  Bg^tip  for  a  fully  magnetized 

nickel  sphere  of  55  /zm  radius.  Gontours  are  cylindrically  symmetric  and  the 
cross-sections  with  the  x-z  plane  are  shown  in  the  figure.  Sweeping  B^xt  to 
higher  field  has  the  effect  of  moving  the  resonant  slice  away  from  the  tip  (to 
contours  of  smaller  Bg^tip  values).  Gontours  of  dBg^tipfdz,  which  determines 
the  thickness  of  the  resonant  slice  and  hence  the  z  resolution  of  MRFM,  are 
also  shown  in  Fig.  5.12.  Contours  of  dBg^npjdx,  which  relates  to  the  lateral 
spatial  resolution  of  MRFM,  are  shown  in  Fig.  5.13. 

As  Bext  is  swept  from  lower  field  with  a  fixed  tip-sample  distance,  initially 
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Bext  (G) 


Figure  5.10:  Force-detected  EPR  spectra  of  a  ~100-/im-thick  phosphorus- 
doped  silicon  sample.  Dopant  concentration  is  ~  8  x  10^®  cm  Arrows 
indicate  a  peak  that  depends  on  the  tip-sample  distance,  (a)  The  tip-sample 
distance  was  ~  80  /im.  (b)  The  tip-sample  distance  was  ~  60  fj,m.  (c)  The 
tip-sample  distance  was  ~  40  /wm. 
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there  is  no  overlap  between  the  sample  and  the  resonant  slice  that  is  too  close 
to  the  tip.  Therefore,  there  is  no  EPR  force  response  until  B^xt  is  swept  up 
such  that  the  resonant  slice  moves  far  enough  from  the  tip  to  reach  the  sample 
surface.  The  small  peaks  observed  in  Fig.  5.10  are  considered  to  correspond 
to  the  response  when  the  resonant  slice  reaches  the  sample  surface.  For 
smaller  tip-sample  distances,  lower  B^xt  is  reqtdred  for  the  resonant  slice  to 
reach  the  sample,  causing  the  small  peak  to  appear  at  lower  field  eis  shown 
in  Fig.  5.10. 

As  Bcxt  is  swept  further  up,  the  overlap  volume  between  the  sample  and 
the  resonant  slice  increases.  However,  the  force  per  unit  volume,  which  is 
proportional  to  \dBx^tip/dz\,  decreases  at  the  same  time,  and  the  total  EPR 
force  changes  rather  gradually.  The  FM  cyclic  saturation  technique  used  in 
the  experiment  gives  a  differential  response  as  explained  in  Section  3.2.2. 
Therefore,  the  gradual  change  beyond  the  point  where  the  resonant  slice 
reaches  the  sample  surface  does  not  result  in  a  distinctive  peak. 

When  Bext  —  4300  G,  the  resonant  slice  is  far  away  from  the  tip  (>  200 
/im),  and  is  thick  because  the  gradient  is  small  far  from  the  tip.  In  other 
words,  a  large  region  of  the  sample  except  the  proximity  to  the  tip  becomes 
on  resonance.  Although  the  force  per  volume  is  quite  small,  the  large  volume 
of  the  resonance  results  in  rather  a  strong  force.  As  B^xt  is  increased  further, 
the  resonant  slice  disappears  rapidly,  causing  a  distinctive  response  in  EPR 
force  spectra  obtained  by  the  FM  cyclic  saturation  technique.  The  peak 
around  4300  G  does  not  strongly  depend  on  the  tip-sample  distance,  since  it 
is  arising  from  the  region  far  away  from  the  tip. 

Force-detected  EPR  spectra  for  a  uniform  thin  sample  can  be  more  rig¬ 
orously  analyzed.  The  EPR  force  from  spins  in  a  volume  dV  is  given  by 

dF  =  ^J^SpMdV,  (5.6) 
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Figure  5.12:  Contours  of  dBzjdz  from  a  fully  magnetized  nickel  sphere  (r  = 
55  /zm).  The  z  axis  denotes  the  distance  from  the  surface  of  the  sphere. 


Figure  5.13:  Contours  of  dBzjdx  from  a  fully  magnetized  nickel  sphere  ( 
55  am).  The  z  axis  denotes  the  distance  from  the  surface  of  the  sphere. 
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where  fis  is  the  Bohr  magneton,  Bq  c::i  4300  G  is  the  resonant  field  for  12-GHz 
microwave  field  J5i,  and  SpM  is  the  response  for  the  FM  cyclic  saturation 
given  by 


Sfm  =  ^ 


Ir 


dBti 


tip 


dz 


(5.6) 


+  [7{Be.t  +  5tip)  -  +  72B?TiT2  , 

Here,  a>{t)  =  Wo+fi  cos{(jJct)  is  the  frequency  of  Bi  modulated  at  the  cantilever 
resonant  frequency  oJc,  and  denotes  the  first-harmonic  Fourier  component 
at  Uc.  The  total  force  can  be  obtained  by  integrating  dF  over  the  sample 
volume.  We  approximate  the  sample  to  be  infinitely  thin  with  the  dopant 
concentration  per  unit  area  crp  =  NpS,  where  5  is  the  sample  thickness. 
Then,  noting  Btip  and  dBtipIdz  are  cylindrically  symmetric,  we  obtain  the 
total  force 

F  =  r““  Sfm  •  27rx  dx,  (5.7) 

/cb-Z  Jo 

where  Xmax  is  a  reasonably  large  cutoff  length. 

Figures  5.14,  5.15,  and  5.16  show  the  calculation  results  of 

/  Sfm  •  27rx  dx  (5-8) 

Jo 

as  a  function  of  B^xt  for  the  tip-sample  distances  of  150  /im,  100  /xm,  and 
50  /xm,  respectively.  Values  of  7  =  1.76  x  10^  rad  s“^G“^,  Ti  =  T2  =  10“® 
s,  Bi=0.5  G,  u;o/27r=11.84  GHz,  Q/2Tr  =  10  MHz,  x^ax  =  1  mm,  and  Btip, 
dBtipIdz  for  a  fully  magnetized  nickel  sphere  of  55  /Lxm  radius  were  used. 
The  results  show  a  bipolar  peak  relatively  independent  of  the  tip-sample 
distance,  and  a  small  peak  that  moves  towards  lower  field  as  the  tip-sample 
distance  is  decreased,  in  good  qualitative  agreement  with  the  experiments. 
The  peak-to-peak  value  is  on  the  order  of  1,  giving 

(TdIJ,%Bo 


F  ciz 


X  1 


keT 

~  4  X  10“^  dyne 
~  4  X  10-^2  N, 


(5.9) 
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Figure  5.14:  The  simulation  result  for  EPR  force  for  a  thin  sample  and  a 
55-/im-radius  nickel  sphere  tip.  Tip-sample  distance  is  150  /um. 


which  is  about  40  times  larger  than  the  experimental  results.  The  discrepancy 
is  presumably  due  to  the  imcertainty  in  the  parameter  values  and  the  rather 
crude  approximation  that  the  sample  is  infinitely  thin. 

Small  features  appearing  at  the  high-held  side  of  the  bipolar  peak  in 
Fig.  5.10  are  probably  due  to  the  hyperhne  splitting,  although  the  size  and 
phase  of  the  signal  is  not  understood. 
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Figure  5.15:  The  simulation  result  for  EPR  force  for  a  thin  sample  and  a 
55-/xm-radius  nickel  sphere  tip.  Tip-sample  distance  is  100  fxm. 
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Figure  5.16:  The  simulation  result  for  EPR  force  for  a  thin  sample  and  a 
55-/Ltm-radius  nickel  sphere  tip.  Tip-sample  distance  is  50  /zm. 
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5.3.2  High-dopant-concentration  silicon 

Results  for  the  high-concentration  (N^  ~  4  x  10^*  cm“^)  sample  showed  a 
remarkable  difference  from  the  low-concentration  sample  results.  Figure  5.17 
shows  force-detected  EPR  spectra  of  the  high-concentration  sample  for  the 
tip-sample  distances  of  100,  80,  60,  and  40  There  is  no  apparent  de¬ 
pendence  on  the  tip-sample  distance.  Besides  the  bipolar  peak  around  4300 
G,  no  peak  is  observed  downfield  imlike  the  case  for  the  low-concentration 
sample.  The  close-up  of  the  bipolar  peaJc  shows  some  dependence  on  the 
tip-sample  distance,  including  the  peak  height  difference,  the  shift  of  the 
peak,  and  the  relatively  long  tail  towards  lower  field  for  smaller  tip-sample 
distances. 

The  lack  of  the  major  dependence  on  the  tip-sample  distance  is  a  con¬ 
sequence  of  the  nature  of  conduction  electrons.  Unlike  the  localized  elec¬ 
trons  in  the  low-concentration  sample,  the  conduction  electrons  in  the  high- 
concentration  sample  feels  an  “average”  field  from  the  tip.  Since  the  Hamil¬ 
tonian  H,  which  has  the  spatially  dependent  part  arising  from  Btip,  appears 
in  the  form 

j  V’*(r)ff(r)V’(r)  dr  (5.10) 

to  give  the  energies  of  the  electron,  the  effect  of  Btip  is  averaged  over  the 
extended  region  covered  by  the  conduction  electron  wavefunction  ip. 


5.4  EPR  force  imaging  of  DPPH  sample 

To  further  investigate  imaging  characteristics  of  paramagnetic  samples,  EPR 
force  maps  were  obtained  by  scanning  a  non-uniform  sample  in  x-y.  Fig¬ 
ure  5.18(a)  shows  an  optical  micrograph  of  the  sample  consisting  of  two 
diphenylpicrylhydrazil  (DPPH)  particles,  ~5  /um  in  size. 


CHAPTER  5.  EPR  AND  FMR  IMAGING 


84 


B(G) 


Figure  5.17:  Force-detected  EPR  spectra  for  a  high-concentration  (4  x  10^® 
cm“^)  phosphorus-doped  thin  silicon  sample.  Tip-sample  distances  were  ap¬ 
proximately  100,  80,  60,  and  40  /xm. 
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Figure  5.18(b)  shows  a  force-detected  EPR  signal  presented  as  a  gray¬ 
scale  x-y  image  for  B^xt  =  3895  G.  As  explained  in  Ref.  [16],  the  EPR  force 
map  has  a  ring-like  character,  reflecting  the  intersection  of  the  scan-plane 
with  the  pajaboloidal  resonant  slice  created  by  the  tip  field  !/>  z)  -1- 

Bext  —  4300  G).  Two  ring-like  responses  are  observed  in  Fig.  5.18(b),  aris¬ 
ing  from  two  sample  particles.  The  center  of  each  ring  corresponds  to  the 
lateral  position  of  each  particle.  The  diameters  of  the  rings  are  indicative  of 
the  z  position  of  the  particles.  Since  the  ring-like  response  corresponds  to 
the  intersection  of  the  sample  scan-plane  and  the  resonant  slice,  which  is  a 
contour  of  Btip  ~  400  G  as  shown  in  Fig.  5.11,  the  diameters,  ~38  /im  and 
'^68  ^m,  respectively,  indicates  the  tip-sample  distances  of  '^56  jtim  and  51 
//m,  respectively.  The  larger  diameter  means  the  closer  tip-sample  distance. 

Each  ring-like  response  in  Fig.  5.18(b)  contains  the  information  of  the  size 
and  shape  of  each  sample  particle  contributing  to  the  response.  This  fact  is 
more  clearly  demonstrated  using  the  sample  consisting  of  larger  particles, 
shown  in  Fig.  5.19. 

Figure  5.20  shows  EPR  force  maps  for  four  B^xt  values  obtained  by  scan¬ 
ning  the  sample  shown  in  Fig.  5.19.  As  in  the  EPR  force  map  for  the  smaller 
sample  (Fig.  5.18),  the  force  maps  have  a  ring-like  character,  reflecting  the 
intersection  of  the  scan-plane  with  the  paraboloidal  resonant  slice.  Start¬ 
ing  from  a  point-like  response  in  Fig.  5.20(a),  the  radius  of  the  ring-like 
response  increases  as  B^xt  is  increased.  Here,  increasing  B^xt  moves  the  reso¬ 
nant  slice  away  from  the  tip  (to  contours  of  smaller  Btip  values),  and  therefore 
is  roughly  equivalent  to  reducing  the  tip-sample  distance  with  a  conversion 
factor  |(5Btip/5z)“^|.  The  bipolar  nature  of  the  spectral  lineshape  is  mani¬ 
fested  in  the  EPR  force  map  as  the  dark  inner  side  of  the  ring-like  response 
with  a  brighter  outer  side.  It  is  also  noted  that  the  EPR  force  map  is  not 
a  simple  narrow  ring  but  rather  a  broadened  superposition  of  more  complex 
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Figure  5.18:  (a)  Optical  micrograph  of  DPPH  sample,  (b)  EPR  force  map 
obtained  by  scanning  the  DPPH  sample  in  x-y  at  Bext  =  3895  G. 
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Figure  5.19:  Optical  micrograpli  of  DPPH  sample  consisting  of  two  particles. 


88 


CHAPTER  5.  EPR  AND  FMR  IMAGING 


features,  since  the  signal  contains  contributions  from  spins  distributed  across 
the  spatially  extended  sample. 

Line-broadening  due  to  the  contributions  from  the  extended  sample  is 
clearly  seen  in  the  EPR  spectrum  shown  in  Fig.  5.21.  The  spectrum  was 
obtained  by  sweeping  with  a  fixed  sample  position,  near  the  center  of 
the  ring-like  response  in  the  force  map  in  Fig.  5.20.  The  spectrum  has  a 
bipolar  lineshape  as  a  result  of  the  FM  cyclic  saturation  technique,  and  a 
broadened  linewidth,  ~100  G  peak-to-peak. 

To  acquire  a  real-space  image,  a  three-dimensional  EPR  force  map  F{x,  y,  5ext) 
was  measured  on  an  x-y-B^xt  lattice  of  64  x  64  x  33  points  with  an  acquisition 
rate  of  0.9  s  per  data  point  (0.6  s  for  cantilever  frequency  and  interferometer 
spacing  adjustment  and  0.3  s  for  signal  integration).  The  measured  force 
map  can  be  written  using  a  real-space  magnetization  map  M{x,y,  z)  as 

F[x,y,  Bfxt)  ~  JJJ  iV' "I  iV  dy  dz  ,  (5.11) 


where 

Sfm 


l^BlTxT^ 


+  {7[-Bext  +  Bt\p{x'  -x,y'-  y,  z')]  -  u}{t)yT^  +  'y^BlTiT2^ 


X 


dBtipjx'  -  x,y'  -  y,  z') 
dz 


(5.12) 


is  the  response  for  the  FM  cyclic  saturation  as  in  Eq,  (5.6).  We  approximate 


that 


Btip(a:,y,z)  ~  Btip(s,y,zo)  +  -  ^o),  (5.13) 


where  a  is  a  fixed  representative  value  of  dBupIdz.  Then  the  B^xt  scan  is 
equivalent  to  the  z  scan  with  the  conversion  B^xt  =  —cxz,  and  Eq.  (5.11)  can 
be  written  as  a  convolution  integral 


Fequiv(®j  y>  ^)  — 


F{x,y,-az)  (5.14) 

jjj  M{x', y',  z')h{x  -x\y-  y',  z  -  z')  dx'  dy'  dz', 
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Figure  5.20:  EPR  force  maps  obtained  by  scanning  the  DPPH  sample  in  x-y 
at  fixed  values  of 
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Figure  5.21:  Force-detected  EPR  spectrum  for  the  DPPH  sample. 
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where 


h{x  -  x\y  -y\z-  z')  = 

j:( _ _ 

VI  +  {7[Btip(a;'  -  x,y' -  y,  zq)  +  a{z'  -  z)  -  o-zq]  - 
dBti^{x'  -  x,y'  -  y,ZQ) 
dzQ 


We  can  reconstruct  the  real-space  M{x,y,z)  using  the  Wiener  filtering  de- 
convolution  technique  described  previously  by  Ziiger  et  af.[15,  16,  17]  such 


M  =  , 

\h\^  +  C 


(5.16) 


where  M,  Fequiv,  and  h  are  the  complex-valued  three-dimensional  Fourier 
transforms  of  M,  Fequiv,  and  h,  respectively,  and  (7  is  a  constant  that  is 
adjusted  according  to  the  signal-to-noise  ratio  of  the  data[43].  The  real-space 
M  is  then  fotmd  by  the  inverse  Fourier  transform  of  M. 

The  point-response  function  h  for  the  deconvolution  was  obtained  using 
calculated  values  of  the  field  from  a  fully  magnetized  nickel  sphere  of  55  fxm 
radius  at  zq  =  73  (xm  as  a  representative  position  and  a  theoretical  formula 
for  EPR  response  based  on  the  estimated  values  of  Bi,  spin  relaxation  times 
for  DPPH,  and  the  FM  deviation  of  Bi  used  in  the  experiment. 

A  constant  density  surface  of  the  reconstructed  real-space  three-dimensional 
magnetization  is  shown  in  Fig.  5.22,®  showing  two  main  particles  consistent 
with  the  optical  micrograph  in  Fig.  5.3.  Evidence  of  the  difference  in  z 
position  of  the  two  particles  can  be  seen  in  the  force  map,  especially  in 
Pig.  5.20(c),  as  the  difference  of  diameters  of  ring-like  responses  arising  from 
the  two  main  particles.  The  spatial  resolution  of  the  reconstructed  image  is 


^Image  processing  software  PV-WAVE  was  used  in  the  deconvolution. 
®IBM  Data  Explorer  was  used  to  render  the  image. 
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estimated  to  be  roughly  5  fim.  Although  overall  shape  and  size  of  the  recon¬ 
structed  objects  appear  to  be  a  good  representation  of  the  sample  particles, 
the  image  seems  to  be  somewhat  distorted  and  the  detailed  structures  are  not 
likely  to  be  real.  These  artifacts  are  due  to  the  finite  signal-to-noise  ratio  and 
the  limitation  of  the  approximation  that  the  Bext  scan  is  equivalent  to  the 
z  scan.  We  would  expect  higher  quality  reconstruction  if  we  had  physically 
scanned  the  tip-sample  distance  rather  than  Bext-  This  was  not  possible  in 
our  present  apparatus  due  to  the  limited  z-scan  range  of  the  piezoelectric 
tube. 


5.5  FMR  force  imaging 

Force  detection  of  ferromagnetic  resonance  (FMR),  recently  demonstrated 
by  Zhang  et  o/.[44],  was  found  to  have  fundamentally  different  imaging  char¬ 
acteristics  compared  to  the  above  paramagnetic  case.  FMR  imaging  was 
studied  for  a  sample  consisting  of  a  small  yttrium  iron  garnet  (YIG)  plate, 
23  fim  X  23  fjLm  and  3  fim  thick.  The  sample  was  prepared  using  pho¬ 
tolithography  and  ion  milling  techniques,  starting  from  a  YIG  thin  film  on  a 
gadolinium  gallium  garnet  (GGG)  substrate.®  The  sample  was  oriented  such 
that  Bext  was  normal  to  the  plate. 

Figure  5.23(a)  shows  a  force-detected  FMR  spectrum  obtained  at  a  fixed 
x-y-z  position  while  B^xt  was  swept.  The  sample  position  was  roughly  50 
/im  below  the  magnetic  tip  and  50  /xm  away  from  the  axis  of  the  tip  lat¬ 
erally.  Despite  the  large  lateral  field  gradient  from  the  tip  at  this  position, 
iVBtipl  Ci;  10  (see  Figs.  5.12  and  5.13),  the  spectrum  shows  no  signifi¬ 

cant  broadening  in  contrast  to  the  force-detected  EPR  spectrum  in  Fig.  5.21. 


®  We  thank  R.  Fontana,  L.  Lane,  and  R.  Campbell  for  the  fabrication  of  the  YIG  sample. 
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top  view 


side  view 


Figure  5.22:  Reconstructed  three-dimensional  magnetization  image.  Con¬ 
stant  density  surfaces  are  depicted  from  two  viewing  orientations. 
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The  peak-to-peak  linewidth  in  Fig.  5.23(a)  is  roughly  8  G,  much  less  than 
one  would  expect  from  the  sample  size  and  the  field  gradient  (23  fim  x  10 
G///m  =  230  G). 

The  lack  of  broadening  is  clearly  shown  in  the  FMR  force  map.  Fig¬ 
ure  5.23(b)  shows  the  FMR  force  signal  obtained  by  scanning  the  sample  in 
x-y  with  fixed  z  and  B^xt  =  4418  G.  The  result  shows  a  ring-like  response  as 
in  the  EPR  force  map;  however,  the  width  of  the  ring  is  narrow  despite  the 
finite  sample  size.  The  difference  is  clear  if  one  compaxes  the  FMR  results 
with  the  EPR  force  maps  in  Fig.  5.20,  where  a  broadened  and  more  complex 
response  was  observed  due  to  the  finite  size  aind  shape  of  the  sample. 

The  narrow  spectral  width  is  a  consequence  of  the  nature  of  FMR.  In 
FMR,  spins  are  exchange  coupled  to  one  another  and  each  spin  does  not 
act  independently.  In  fact,  in  the  uniform  mode,  which  we  believe  we  are 
observing,  all  spins  are  precessing  uniformly.  This  makes  spins  in  the  finite 
size  sample  act  like  a  single  collective  spin  ensemble  and  the  FMR  force  map 
resembles  the  response  from  a  point-like  sample. 

In  addition  to  the  main  spectral  line  in  Fig.  5.23(a),  two  small  magnetic 
resonance  lines  were  observed  downfield  from  the  main  signal.  These  are 
higher-order  magnetostatic  modes  as  are  observed  in  conventional  FMR  for 
normally  magnetized  YIG  disks[45].  In  the  FMR  force  map.  Fig.  5.23(b), 
these  modes  appear  as  faint  concentric  rings  with  larger  diameters  than  the 
main  signal.  The  concentric  rings,  including  the  main  signal,  represent  con¬ 
tours  of  Btip  for  constant  values  that  correspond  to  the  resonant  fields  of  each 
mode  (offset  by  the  uniform  field  jBext)-  Larger  diameter  rings  correspond  to 
contours  of  lower  field  values. 

Measurement  with  lower  microwave  power  (smaller  Bi)  revealed  addi¬ 
tional  spectral  features  within  the  main  spectral  line,  as  shown  in  Fig.  5.24(a). 
We  believe  that  these  are  the  same  structures  as  observed  by  Dillon[45]  us- 
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Figure  5.23:  (a)  Force-detected  FMR  spectrum  for  a  small  YIG  plate,  show¬ 
ing  a  narrow-linewidth  main  signal.  Arrows  indicate  small  spectral  features 
observed  due  to  magnetostatic  modes,  (b)  FMR  force  map  at  B^xt  —  4418 
G,  showing  a  sharp  ring-like  response.  Arrows  indicate  concentric  ring-like 
responses  due  to  magnetostatic  modes. 
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ing  conventional  FMR  and  proposed  to  be  due  to  magnetoelastic  coupling 
of  magnetostatic  and  acoustic  modes.  In  the  FMR  force  map  with  fixed 
Bext  =  4342  G,  shown  in  Fig.  5.24(b),  these  spectral  features  are  not  seen 
as  simple  concentric  ring-like  responses  as  in  the  case  for  the  magnetostatic 
modes.  Rather,  these  features  seem  to  be  spiraling  across  the  main  response 
depending  on  the  x-y  position,  suggesting  that  the  magnetoelastic  interac¬ 
tion  is  not  uniform  over  the  entire  sample.  This  result  illustrates  the  ability 
of  MRFM  to  visualize  complex  magnetic  resonance  phenomena  on  the  mi¬ 
crometer  scale. 


5.6  Magnetic  field  mapping  using  force- detected 
FMR 

The  sharp  spectral  features  of  FMR  can  be  used  to  make  quantitative  mea¬ 
surements  of  magnetic  field  above  a  ferromagnetic  sample.  As  shown  in 
Fig.  5.25,  the  idea  is  to  use  force  detection  to  measure  ferromagnetic  res¬ 
onance  of  a  small  yttrium-iron-garnet  (YIG)  particle  while  the  particle  is 
mounted  on  a  cantilever  and  scanned  over  a  sample.  Figure  5.26,  shows  the 
force-detected  spectrum  from  a  micrometer-size  gallium-doped  YIG  particle. 

A  number  of  resonance  peaks  appear  at  well-defined  field  values.  The  peaks, 
which  are  due  to  vaxious  magnetostatic  spin  wave  modes  in  the  YIG  parti¬ 
cle,  can  be  used  as  field  markers  when  the  YIG  particle  is  scanned  over  a 
sample  that  generates  an  inhomogeneous  field  pattern.  As  a  first  example 
of  quantitative  field  mapping,  the  YIG  particle  was  scanned  over  a  sample 
consisting  of  a  20m  diameter  nickel  sphere.  The  image  in  Fig.  5.27  shows 
a  series  of  bright  rings  centered  over  the  nickel  particle.  Each  ring  corre¬ 
sponds  to  a  specific  contour  of  constant  field  at  a  distance  approximately  20 
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Figure  5.24:  (a)  Force-detected  FMR  spectrum  obtained  at  lower  microwave 
power.  Small  spectral  features  believed  to  be  due  to  magnetoelastic  coupling 
are  observed  within  the  main  line,  (b)  FMR  force  map  at  Bcxt  =  4342  G 
with  the  same  microwave  power.  Note  the  spiral-like  features  superimposed 
on  the  main  circular  response. 
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Microwave  Coil  - 1.7  GHz 


Figure  5.25:  Apparatus  to  quantitatively  map  the  magnetic  field  emanating 
from  a  sample.  A  Ga-doped  YIG  particle  mounted  on  the  cantilever  serves  as 
a  local  probe  of  magnetic  field  emanating  from  a  nearby  sample.  A  cantilever 
vibration  signal  is  generated  whenever  the  field  from  the  sample  matches  the 
field  for  ferromagnetic  resonance  of  the  YIG  particle. 
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Figure  5.26:  Force-detected  ferromagnetic  resonance  spectrum  of  Ga-doped 
YIG  particle. 
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Figure  5.27:  MRFM  image  showing  contours  of  constant  field  above  a  20 
micron  diameter  nickel  sphere. 
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fim  above  the  sample.  The  spatial  resolution  is  limited  only  by  the  size  of 
the  YIG  probe  particle.  By  using  smaller  probe  particles,  extension  to  sub¬ 
micron  resolution  should  be  possible.  In  addition  to  obvious  applications  for 
characterizing  magnetic  recording  heads,  etc.,  the  technique  should  be  use¬ 
ful  for  characterizing  the  field  emanating  from  magnetic  tips  developed  for 
single-spin  MRFM  experiments. 


Chapter  6 

Single  Spin  Experiment 

6.1  Experimental  configuration 

The  ultimate  goal  of  IBM’s  MRFM  research  is  to  achieve  single  spin  sen¬ 
sitivity  with  angstrom  spatial  resolution  in  three-dimensions.  Considerable 
effort  during  the  ONR  contract  h<is  been  devoted  to  developing  the  apparatus 
and  techniques  for  single  spin  detection.  The  emphasis  has  been  on  electron 
spins  rather  than  nuclear  spins  because  of  the  larger  magnetic  moment  of  the 
electron. 

Figure  6.1  shows  a  schematic  diagram  of  IBM’s  single  spin  experiment. 
The  experiment  uses  an  ultrathin  single-crystal  silicon  cantilever  with  a  mag¬ 
netic  tip  positioned  within  10-50  nm  of  the  surface  of  a  polished  Si02  sample. 
The  sample  is  prepared  by  gamma  irradiation  so  that  it  contains  a  known 
density  of  dangling  bond  defects  with  unpaired  electron  spins  (E’  centers) 
throughout  its  volume.  The  cantilever  is  oriented  perpendiculax  to  the  sam¬ 
ple  in  order  to  allow  the  tip  to  approach  the  sample  surface  to  within  10  nm 
without  being  snapped  into  the  surface  by  electrostatic  and  van  der  Waals 
forces.  The  sharp  end  of  the  cantilever  is  coated  along  one  sidewall  with 
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Figure  6.1:  Experiment  to  detect  individual  electron  spins  (E’  centers)  in 
Si02.  The  1  A-thick  resonant  slice  is  shown  by  the  dotted  lines. 
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60nm  of  cobalt  to  form  a  magnetic  tip.  Angstrom-scale  vibration  of  the  can¬ 
tilever  is  detected  by  a  fiberoptic  interferometer.  The  apparatus  is  operated 
at  4K  so  as  to  improve  the  force  resolution  of  the  cantilever  and  to  maximize 
the  spin-lattice  relaxation  time  of  the  spins. 

As  is  usual  for  MRFM  experiments,  the  magnetic  field  from  the  tip  plus 
an  external  field  (in  this  case  from  a  superconducting  magnet)  defines  a 
resonant  slice  wherein  magnetic  resonance  occurs.  Electron  spins  located 
within  the  resonant  slice  undergo  magnetic  resonance  stimulated  by  the  12 
GHz  microwave  magnetic  field  generated  by  a  microstripline  resonator.  The 
microstripline-sample  combination  is  mounted  on  a  piezoelectric  tube  scan¬ 
ner,  allowing  the  tip  to  be  scanned  with  respect  to  the  sample  in  three  di¬ 
mensions. 

The  perpendicular  cantilever  configuration  in  our  experiment  responds  to 
lateral  forces.  The  lateral  force  due  to  the  interaction  of  the  tip  field  gradient 
with  the  target  spin  is  given  by 

F,  =  (6.1) 

where  is  the  magnetic  moment  of  the  electron  (9.3  x  10“^^  emu  or  9.3  x 
10-24  j/T). 

Because  fie  is  so  small,  it  is  important  to  have  a  large  field  gradient. 
For  our  present  tip  fabrication  techniques,  a  gradient  up  to  dBzjdx  =  10 
Gauss/ A  may  be  achievable.  K  this  gradient  is  achieved,  then  the  resulting 
force  is  93  aN  or,  if  the  force  oscillates  sinusoidally,  66  aN-rms.  (laN  = 
1  attonewton  =  10”^®  newton.)  Thus,  despite  the  large  field  gradient,  the 
resulting  magnetic  force  is  very  small  indeed.  For  comparison,  this  force  is 
a  million  times  smaller  than  the  piconewton  forces  detected  by  conventional 
atomic  force  microscopes. 
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6.2  Ultrasensitive  silicon  cantilevers 


Because  of  the  smallness  of  the  force  to  be  detected,  considerable  effort  has 
been  devoted  to  developing  the  capability  to  detect  attonewton  forces.  To 
enhance  our  capability  in  the  area  micromechanical  fabrication,  a  strong 
collaboration  has  been  established  with  Prof.  Tom  Kenny  and  his  student 
Tim  Stowe  at  Stanford  University.  Working  together,  we  have  demonstrated 
the  most  sensitive  force  detection  ever  achieved  with  mechanical  sensors. 

As  discussed  in  Chapter  2,  the  minimum  force  detectable  in  a  bandwidth 
Ai/  is  limited  by  the  thermal  motion  of  the  cantilever  according  to 

The  factor  k/ojc  can  be  written  in  terms  of  cantilever  geometry.  For  a  rect¬ 
angular  cantilever  with  thickness  t,  width  w  and  length  Z,  the  minimum 
detectable  force  is 

(6.3) 

where  E  is  the  Young’s  modulus  and  p  is  the  density  of  the  cantilever  mate¬ 
rial. 

The  above  equation  suggests  the  strategy  we  have  adopted  to  achieve 
ultrasensitive  force  detection:  fabricate  high  Q  cantilevers  that  are  thin, 
narrow,  and  long.  In  particular,  we  have  succeeded  in  fabricating  single¬ 
crystal  silicon  cantilevers  as  thin  as  500A  and  as  long  as  0.4  mm.  Single¬ 
crystal  silicon  is  a  good  choice  as  the  cantilever  material  since  it  has  low 
intrinsic  mechanical  loss  (i.e.,  high  Q)  and  very  low  internal  stress,  allowing 
long,  thin  structures  to  be  fabricated  without  unwanted  curling. 

As  shown  in  Fig.  6.2,  the  cantilever  processing  begins  with  (lOO)-oriented 
silicon-on-insulator  (SOI)  wafers  with  a  top  silicon  layer  thickness  of  ISOOA. 
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Side  view  Blow  i?)  of  the  tip  region  (top  view) 


I  I  Oxide 

Low  stress  alicm  nitride 


Figure  6.2:  Processing  steps  for  fabricating  ultrathin  silicon  cantilevers  with 
integrated  in-plane  ferromagnetic  tips. 
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This  top  silicon  layer  is  first  thinned  down  to  a  desired  thickness  by  growing 
a  thermal  oxide,  which  is  then  removed  using  a  standard  buffered  oxide  etch 
(BOE).  Next,  the  cantilever  and  the  tip  are  defined  using  a  two  step  mask  and 
etch  process  [Fig.  6.2(b)].  The  first  photoresist  mask  defines  the  cantilever 
outline  with  a  tip  rounded  to  0.5  ^m,  the  limit  of  our  optical  lithography. 
This  cantilever  pattern  is  etched  into  the  silicon  using  a  SFe  plasma  etch. 
This  mask  is  stripped  and  a  second  photoresist  mask  containing  the  same 
pattern  is  applied  again,  but  shifted  to  the  right  of  the  first  pattern  by  1 
lj,m.  A  tip  forms  at  the  intersection  of  the  right  sidewall  formed  during  the 
first  etch  and  the  left  sidewall  of  the  second  etch,  yielding  a  tip  radius  below 
500-A. 

Once  the  cantilevers  have  been  patterned,  a  7500-A-thick  layer  of  low 
temperature  oxide  (LTO)  and  a  l-/fm-thick  layer  of  low  stress  nitride  are  de¬ 
posited  for  topside  protection.  Openings  in  the  backside  LTO  and  nitride  are 
patterned  and  the  exposed  silicon  is  etched  anisotropically  in  20%  tetram- 
ethylammonium  hydroxide  (TMAH).  After  the  backside  etch  is  complete  [see 
Fig.  6.2(c)],  the  top  nitride  layer  is  removed  using  a  SFg  plasma  etch  and  the 
oxide  membrane  encapsulating  the  cantilevers  is  removed  using  BOE. 

The  wafer  is  transferred  under  liquid  from  BOE  to  water  and  then  from 
water  to  methanol,  followed  by  critical  point  drying  in  high  pressure  liquid 
CO2.  This  CO2  drying  step  avoids  meniscus  forces  that  can  easily  breaJr 
or  permanently  curl  the  cantilevers.  Finally,  silicon  dies  defined  during  the 
backside  etch  are  cleaved  [Fig.  6.2(d)].  Each  die  contains  cantilevers  that  are 
4-6  ixm  wide  and  20-400  fxm  long.  Yields  over  80%  have  been  achieved  for 
cantilever  arrays  as  thin  as  500A. 

Figure  6.3  shows  a  finished  220-/um-long,  600-A-thick  cantilever.  The 
cantilever  design  includes  a  20-/ixm-wide  paddle  for  easy  alignment  of  the 
fiberoptic  interferometer  to  the  cantilever.  Tip  radii  were  evaluated  by  elec- 
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Figure  6.3:  A  600  A-thick,  220  ^m-long  single-crystal  silicon  cantilever  with 
in-plane  tip. 

tron  microscopy  and  found  to  be  routinely  in  the  30—60  nm  range. 


6.3  Attonewton  force  detection 

To  evaluate  the  force  resolution  of  a  600  A  thick,  220-//m-long  cantilever, 
its  motion  was  measured  in  vacuum  at  both  room  temperature  and  4.8K 
using  a  fiber  optic  interferometer.  The  thermal  motion  of  the  cantilever  at 
room  temperature  was  measured  to  be  250  A-rms  which,  from  the  equiparti- 
tion  theorem  (equation  (2.1)),  implies  a  spring  constant  of  6.5  x  10"®  N/m. 
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This  estimated  spring  constant  is  within  10%  of  the  value  predicted  by  finite 
element  analysis. 

At  4.8K,  the  cantilever  thermal  vibration  weus  reduced  to  43A-rms,  some¬ 
what  larger  than  the  theoretical  value  of  32  A-rms  expected  from  T^^^  scal¬ 
ing  of  the  room  temperature  noise.  This  small  discrepancy  may  be  due 
to  laser  heating  of  the  cantilever  or  environmental  vibrations.  Combining 
the  measured  cantilever  vibration  noise  with  the  other  cantilever  parameters 
[k  =  6.5^N/m,  /o  =  1.7  kHz  and  Q  =  6700)  gives  a  theoretical  noise  of  5.6 
aN/Hz^/2 

To  demonstrate  attonewton  force  detection,  an  electrode  located  ~  1  mm 
from  the  cantilever  was  used  to  generate  a  small  electrostatic  force  oscillat¬ 
ing  at  the  cantilever  resonance  frequency.  Figure  6.4  shows  a  time  trace  of 
the  cantilever  vibration  amplitude  as  detected  by  the  interferometer  and  a 
lock-in  amplifier.  When  the  oscillating  electrostatic  force  was  turned  on,  the 
cantilever  reached  a  steady-state  amplitude  of  370-A-rms.  Based  on  the  pre¬ 
viously  determined  values  of  k  and  Q,  this  corresponds  to  a  36  aN-rms  force. 
Using  a  detection  bandwidth  that  was  dominated  by  the  natural  bandwidth 
of  the  cantilever,  a  noise  level  of  3.6  aN  wa.s  measured,  yielding  a  signal-to- 
noise  ratio  of  10. 

The  demonstrated  force  resolution  of  5.6  aN /Hz^/^  is,  to  our  knowledge, 
the  best  ever  achieved  using  a  a  mechanical  sensor.  Even  better  resolution 
could  have  been  achieved  if  higher  Q  (lower  dissipation)  had  been  achieved. 
Studies  of  energy  dissipation  in  thin  microcantilevers  are  ongoing. 


6.4  Magnetic  tips 

The  tip  of  the  silicon  cantilever  is  made  magnetic  by  evaporating  a  thin  cobalt 
film  onto  the  cantilever  sidewall,  as  shown  in  Fig.  6.2(e).  By  evaporating  the 
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Figure  6.4:  Time  trace  of  the  cantilever  vibration  amplitude  showing  noise 
and  response  to  an  oscillating  36  aN-rms  electrostatic  force.  The  cantilever 
vibration  amplitude  is  shown  on  the  right  hand  scale  and  the  corresponding 
force  on  the  left  hand  scale.  A  force  noise  of  3.6  aN  was  measured.  The 
bandwidth  of  the  measurement  was  set  by  the  natural  bandwidth  of  the 
cantilever  (approximately  0.4  Hz). 
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cobalt  preferentially  onto  the  sidewall  of  the  tip,  a  long,  thin  magnetic  “wire” 
is  formed  that  exhibits  strong  shape  anisotropy.  Side-wall  deposition  has  an 
additional  benefit:  it  avoids  problems  with  stress  in  the  deposited  film  which 
would  cause  the  the  cantilever  to  curl  if  the  film  were  deposited  onto  the 
face  of  the  cantilever.  For  most  of  our  tips,  knife-edge  shadow  masking  was 
used  to  confine  the  cobalt  to  the  tip  region,  thus  minimizing  mechanical 
dissipation  effects  due  to  the  presence  of  the  cobalt. 

Magnetic  characteristics  of  the  tips  have  been  investigated  in  several  ways. 
First,  the  tips  were  used  for  magnetic  force  microscopy  (MFM)  of  bit  transi¬ 
tions  on  magnetic  recording  disks.  High  quality  MFM  images  were  obtained. 
In  addition,  information  about  the  total  moment  and  anisotropy  of  the  tip 
was  obtained  by  measuring  the  cantilever  frequency  as  function  of  applied 
magnetic  field. 

Figure  6.5  shows  a  hysteresis  loop  obtained  by  measuring  the  natural 
frequency  of  the  cantilever  as  a  function  of  the  magnetic  field  applied  along 
the  longitudinal  axis  of  the  cantilever.  If  we  start  at  high  field,  where  the  tip 
magnetization  is  saturated,  and  follow  the  frequency  as  the  field  is  decreased 
towaxd  zero,  we  see  a  smooth  change  in  frequency.  As  the  field  goes  negative, 
there  are  jumps  in  frequency  occurring  at  approximately  -900  Gauss  and  - 
6000  Gauss.  These  frequency  jumps  represent  sudden  switching  of  the  tip 
magnetization.  The  large  jump  at  -6000  Gauss  is  believed  to  be  due  to 
magnetic  switching  of  the  cobalt  located  along  the  sidewall  of  the  cantilever, 
where  the  bulk  of  the  cobalt  is  located.  This  magnetic  material  exhibits  a 
large  coercivity  because  of  the  strong  shape  anisotropy.  The  small  jump  at  - 
900  Gauss  is  believed  to  be  due  to  the  small  amount  of  material  inadvertantly 
deposited  onto  the  face  of  the  cantilever.  The  coercivity  is  much  smaller 
because  there  is  very  little  shape  anisotropy. 

The  portion  of  the  hysteresis  loop  from  15000  Gauss  to  -900  Gauss  is 
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Hysteresis  Loop  of  Cobalt-coated  Tip 

Cantilever  frequency  magnetometry 


Figure  6.5:  Hysteresis  loop  of  cantilever  frequency  vs.  applied  field. 
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reversible  and  shows  very  little  hysteresis  if  the  direction  of  the  field  sweep 
is  reversed.  This  is  an  indication  that  the  tip  magnetization  remains  fully 
saturated  over  this  field  range. 

The  behavior  of  the  reversible  portion  of  the  hysteresis  loop  (+15000 
Gauss  to  -900  Gauss)  can  be  understood  using  a  simple  model  that  combines 
cantilever  motion  with  the  Stoner- Wohlforth  model  of  single  domain  mag¬ 
netic  particles.  Referring  to  Fig.  6.6,  the  model  assumes  that  the  magnetic 
tip  can  be  modeled  as  an  elongated  single  domain  particle  with  saturation 
magnetization  M*  and  uniaxial  anisotropy  in  the  presence  of  a  z-directed 
applied  field.  The  cantilever  is  also  oriented  in  the  z  direction,  but  its  free- 
end  is  allowed  to  bend  as  it  vibrates.  The  angle  a  denotes  the  tilt  angle 
of  the  tip  due  to  the  bending  of  the  cantilever,  a  is  related  to  the  lateral 
displacement  of  the  tip  according  to 

x  =  Lsma,  (6.4) 

where  X  is  an  effective  cantilever  length  that  is  determined  by  the  shape  of 
the  cantilever  vibrational  mode.  The  direction  of  the  longitudinal  axis  of  the 
particle  is  assumed  to  follow  the  direction  of  the  free-end  of  the  cantilever. 
The  angle  of  the  tip  magnetization  with  respect  to  the  longitudinal  axis  of 
the  particle  is  denoted  by  $. 

The  potential  energy  of  the  cantilever-magnetic  particle  combination  can 
be  written  as  the  sum  of  anisotropy,  Zeeman  and  elastic  energy  terms: 

E  =  Ku  sin^  0  —  HMs  cos(q:  —  ^)  +  sin^  a.  (6.5) 

The  angle  Q  of  the  particle  magnetization  is  assumed  to  lie  in  the  direction 
that  minimizes  the  energy.  Accordingly,  we  set  dEfdB  =  0,  yielding 


2Ku  sin  6  cos  9  —  HMs  sin(a  —  0)  =  0. 


(6.6) 


114 


CHAPTER  6.  SINGLE  SPIN  EXPERIMENT 


z  i 


Figure  6.6:  Model  of  single  domain  particle  on  a  vertical  cantilever.  The 
easy  (long)  cixis  of  the  particle  is  assumed  to  be  aligned  with  the  free  end 
of  the  cantilever  which,  during  vibration,  tilts  by  an  angle  a.  The  particle 
magnetization  is  canted  by  an  angle  0  with  respect  to  the  particle  easy  axis 
because  of  a  z-directed  applied  field. 
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Solving  this  for  small  angles  gives 

e=(  y. 

\2K„  +  HMJ 

Defining  the  anisotropy  field  as  Hk  =  2KulMs  results  in 

.  f  H  \ 


\Hk  +  HJ  ^  ' 

Substituting  (6.4)  and  (6.8)  into  (6.5)  and  using  the  small  angle  approxima¬ 
tion  gives  a  spring-like  potential  energy  with  a  modified  spring  constant: 

E  =  i(fc  4- Afc)^c^  (6.9) 

where  the  spring  constant  change  due  to  the  presence  of  the  applied  field  is 

This  simple  model  gives  a  virtually  perfect  fit  with  experimental  data,  as 
can  be  seen  in  Fig.  6.7.  In  this  figure  we  show  as  points  the  measured  change 
in  spring  constant  as  determined  from  the  frequency  data  according  to 


AA:  =  fc  1 


(6.11) 


where  /o  is  the  cantilever  frequency  at  zero  field.  The  solid  curve  in  Fig. 
6.7  represents  a  two  parameter  fit  of  equation  (6.10)  to  the  experimental 
data.  The  fitting  process  uses  the  known  spring  constant  and  length  of  the 
cantilever  and  yields  an  estimate  of  the  moment  of  the  tip  (1.8  x  10“®  emu) 
and  the  anisotropy  field  (8940  Oe). 

The  fact  that  the  tip  exhibits  such  strong  shape  anisotropy  and  remains 


in  a  saturated  magnetic  state  down  to  zero  external  field  is  very  good  news. 
This  should  allow  the  tip  to  produce  the  maximum  field  and  field  gradient 
possible  for  the  given  geometry.  To  further  evaluate  the  details  of  the  tip  field 
near  the  end  of  the  tip,  we  have  recently  given  samples  of  our  tips  to  Prof. 


M.  Scheinfein  at  Arizona  State  University  for  characterization  by  electron 
holography. 
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Figure  6.7:  The  points  show  the  increase  in  cantilever  stiffness  measured  as  a 
function  of  applied  field  for  the  reversible  portion  of  the  hysteresis  loop.  The 
solid  curve  is  a  plot  of  equation  (6.10)  using  the  fitting  parameters  shown. 
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6.5  Single  spin  experiments 

The  single  spin  detection  apparatus  is  essentially  complete  and  has  been  used 
in  several  attempts  to  detect  single  spin  magnetic  resonance.  Although  no 
single  spin  signal  has  yet  been  seen,  we  have  gained  valuable  experience  in 
operating  our  ultrasensitive  cantilevers  in  close  vicinity  to  sample  surfaces. 

Several  different  methods  for  detecting  the  single  spin  signal  are  under 
investigation.  Here  are  the  three  leading  contenders; 

Method  1  -  Cyclic  adiabatic  inversion  induced  by  frequency  mod¬ 
ulation  of  the  microwave  field 

A  frequency  modulated  microwave  field  is  used  to  induce  cyclic  adiabatic 
inversion  of  a  resonant  spin.  The  cyclic  inversion  of  the  spin  generates  an 
oscillating  force  on  the  cantilever,  causing  it  to  vibrate. 

Advantage:  This  is  a  “zero-baseline”  experiment,  which  readily  allows 
the  theoretical  thermal  noise  limit  to  achieved. 

Disadvantages:  a)  It  requires  a  heavily  damped  cantilever  so  that  can¬ 
tilever  vibrations  due  to  thermal  noise  and/or  the  spin  signal  don’t  exceed 
0.1  A  .  This  is  required  so  that  the  variation  in  the  tip  field  due  to  the  tip 
vibration  As  and  the  tip  field  gradient  doesn’t  disrupt  the  magnetic  reso¬ 
nance  of  the  spin,  (i.e.,  As  dBzjdx  <C  D/7,  where  D  is  the  peah  frequency 
deviation  of  the  frequency  modulation.)  b)  The  cw  microwave  field  can  lead 
to  heating  of  the  experiment,  c)  Spin  relaxation  time  for  spins  constantly 
undergoing  cyclic  adiabatic  inversion  may  be  short  (~  40  ms). 

Method  2  -  Tip-induced  cyclic  adiabatic  inversion 

Cyclic  adiabatic  inversion  is  induced  by  using  a  cw  microwave  field  of  constant 
frequency  and  vibrating  the  cantilever.  The  cyclic  inversion  of  the  target 
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spin  induced  by  the  oscillating  tip  field  will  result  in  a  small  shift  in  the 
eigenfrequency  of  the  cantilever. 

Advantage:  Effectively  enlarges  the  thickness  of  the  resonant  slice  and 
avoids  problem  of  keeping  the  target  spin  resonant  in  the  presence  of  tip 
thermal  vibrations. 

Disadvantages:  a)  Detection  of  spin-induced  cantilever  frequency  shift 
must  be  distinguished  from  much  larger  frequency  shifts  induced  by  mis¬ 
cellaneous  tip-surface  interactions,  b)  The  cw  microwave  field  can  lead  to 
heating  of  the  experiment,  c)  Spin  relaxation  time  for  spins  constantly  un¬ 
dergoing  cyclic  adiabatic  inversion  may  be  short  (~  40  ms). 

Method  3  -  Sub-harmonic  tip-induced  adiabatic  inversion^ 

The  microwave  field  is  off  most  of  the  time  but  is  periodically  pulsed  on  for 
half  of  a  cantilever  vibration  cycle  (e.g.  once  per  hundred  cantilever  cycles). 
The  combination  of  the  pulsed  microwave  field  with  the  field  sweep  caused  by 
the  cantilever  motion  inverts  the  spin.  Interaction  between  the  spin  and  the 
magnetic  tip  shifts  the  frequency  of  the  cantilever,  which  alternates  between 
a  positive  and  negative  frequency  shift  each  time  the  microwave  pulse  inverts 
the  spin. 

Advantage:  Low  duty  cycle  of  the  microwave  pulse  increases  spin  relax¬ 
ation  time  and  reduces  heating  of  the  experiment. 

Disadvantage:  Frequency  shift  is  smaller  than  in  Methods  1  and  2. 

We  are  presently  exploring  all  three  methods  to  better  determine  their 
relative  merits.  Method  3  appears  to  be  the  front-runner  at  the  present  time. 


^This  method  was  suggested  to  us  by  J.  Sidles 
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6.6  Key  outstanding  issues 

Several  key  issues  in  the  single  spin  experiment  remain  to  be  better  under¬ 
stood.  These  issues  axe: 

1.  Is  the  gradient  from  our  tip  large  enough?  If  we  don’t  achieve  at  least 
5  Gauss/ A,  we  may  not  have  sufficient  force  resolution  to  detect  the 
signal.  We  are  presently  trying  to  evaluate  this  gradient  using  electron 
holography. 

2.  What  is  the  lifetime  of  the  spin  quantum  state  during  measurement? 
One  concern  is  anomalously  rapid  relaxation  induced  by  the  proximity 
of  the  ferromagnetic  tip.  This  is  difficult  to  evaluate,  but  should  be 
improved  by  using  tips  with  large  uniaxial  einisotropy.  If  the  relaxation 
rate  is  too  fast  (e.g.,  >  10  sec“^),  then  the  detection  bandwidth  needs 
to  be  commensurately  wide,  resulting  in  increased  noise. 

3.  Are  E'  centers  near  the  surface  of  our  Si02  sample  similar  in  nature  to 
bulk  E'  centers  in  terms  of  number  density,  relaxation  times,  etc.  We 
know  a  lot  about  E'  centers  in  the  bulk,  but  not  near  the  surface. 

4.  The  tip-sample  interaction  due  to  van  der  Waals  and  electrostatic  forces 
is  much  stronger  than  the  single  spin  magnetic  force.  Will  the  single 
spin  signal  be  overwhelmed  by  these  other  forces? 

Resolving  these  remaining  issues  will  be  a  key  part  of  the  follow-on  ONR 


contract. 
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